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A DATA REDUCTION TECHNIQUE AND ASSOCIATED COMPUTER 
PROGRAM FOR OBTAINING VEHICLE ATTITUDES WITH 
A SINGLE ONBOARD CAMERA 

By Richard J. Bendura and Philip G. Renfroe* 

Langley Research Center 

SUMMARY 

A detailed discussion of the application of a previously developed method to deter- 
mine vehicle flight attitude using a single camera onboard the vehicle is presented with 
emphasis on the digital computer program format and data reduction techniques. Appli- 
cation requirements include film and Earth- related coordinates of at least two landmarks 
(or features), location of the flight vehicle with respect to the Earth, and camera charac- 
teristics. Included in this report are a detailed discussion of the program input and out- 
put format, a computer program listing, a discussion of modifications made to the initial 
method, a step-by-step basic data reduction procedure, and several example applications. 
The computer program is written in FORTRAN IV language for the Control Data 6000 
series digital computer. 

INTRODUCTION 

A postflight photogrammetric method was previously devised for determining a con- 
tinuous history of vehicle flight attitudes (Euler angles) using only film data from a single 
onboard camera along with a ground track of the vehicle. A discussion of method 
requirements, assumptions, mathematical relationships, and results is presented in ref- 
erence 1. The method is based on work presented in references 2 and 3. Results from 
other applications and comparisons with statistical trajectory reconstruction techniques 
are included in references 4 to 6. However, detailed discussions of method application 
or data reduction system capabilities are not included in the reference documents. 
Numerous inquiries have been received concerning method application, particularly from 
the standpoint of computer programing. In addition, several modifications to the program 
which make the method more versatile and streamline its use have been made for appli- 
cation to the Viking Balloon Launched Decelerator Test (BLDT) program data reduction 
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(refs. 4 and 5). The purpose of this paper is to present details concerning method appli- 
cation, including an updated program listing and example applications, and to discuss 
current capabilities and restrictions of the associated Langley Research Center data 
reduction system. This paper is intended as a user's guide for future applications of 
the method. 

SYMBOLS 

Measurements are presented in both SI units and U.S. Customary Units. The 
measurements and calculations were made in U.S. Customary Units. 

c focal length of camera 

h height above sea level of camera, meters (feet) 

I point in image space 

K constant of proportionality between object space coordinate system and 

image space coordinate system 

n number of observations 

Q point in object space 

R geocentric earth radius, meters (feet) 

X,Y, Z object space coordinate axes with origin at focal point; also, distance along 

these axes, meters (feet) 

X,Y, Z distance from origin of geocentric coordinate axes, meters (feet) 

X E ,Y E ,Z E Earth-fixed axis system 

Xj,Yj, Zj image space coordinate axes with origin at focal point 

x,y measured image coordinates in image space coordinate system 

x ,y coordinates of point of intersection between focal plane and focal axes in 

Jr r 

image space coordinate system; center of frame 
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a, 1 3 


angles of attack in pitch and yaw, respectively, degrees 


Act 


change in attitude angle, degrees 


error in x 


error in y 


V 

A 

x 'l> M l> ' 

x 2> V-2’ v 2 
x 3’^3> v 3 

a l> a 2’ a 3 

<f> 


total angle of attack, degrees 


longitude of camera measured positive east from Greenwich, degrees 


directional cosines of image space coordinate system relative to object 
space coordinate system 


camera azimuth, pitch, and roll Euler angles relative to Earth-fixed axes, 
degrees 

geodetic latitude of camera measured positive north, degrees 

geocentric latitude of camera measured positive north, degrees 

vehicle yaw, pitch, and roll Euler angles relative to Earth-fixed axes, 
degrees 


Subscripts: 


Superscript: 


for camera i = 1 and for object space points i = 2, . . ., n + 1 


initial conditions 


vehicle or spacecraft 


denotes transpose of matrix 
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CAMERA- VEHICLE ORIENTATION METHOD 


In order to determine the orientation of a flight vehicle (or spacecraft) with respect 
to the Earth by using this photogrammetric method, the following information is required: 
location of at least two landmarks on each frame of the film, Earth- related coordinates 
of these landmarks, location of the flight vehicle with respect to Earth (as determined by 
radar, for example), orientation of the camera within the vehicle, camera focal length, 
and lens distortion characteristics. Atmospheric refraction corrections are not included 
because these corrections would have been insignificant for previous applications. 

Digital computer techniques can then be applied to determine the relationship between 
the two coordinate systems defined by the film frame (image space) and the Earth (object 
space). For the convenience of the potential user, the mathematics of the relationship 
(taken directly from ref. 1) are included as appendix A. 

Modifications 

For the previous applications of the method, two separate digital computer pro- 
grams were employed to operate on the raw data (image space identification points read 
from the film) before the final camera or vehicle orientation angles (Euler angles, figs. 1 
and 2) were determined. Also a third program was required to operate on the vehicle 
Euler angle (i//,0,$) data and produce vehicle angles (a,fi,ri) relative to the wind. In the 
modified version, presented in this paper, the three programs have been combined into 
one without any reduction in data output and with significant improvement in user and 
computer efficiency. The new version also includes simplified input-output procedures 
and a multiple job processing capability. Another significant modification, to be dis- 
cussed in the following section, involves converting unknown surface features or non- 
permanent objects appearing on the film into usable "landmarks." 

Method Application 

Application of the method is discussed in terms of procedures by use of the basic 
data reduction and computer systems currently available at the NASA Langley Research 
Center and as used to determine test vehicle motions from the Viking Balloon Launched 
Decelerator Test (BLDT) Program (refs. 4 and 5). 

Landmark selection .- The initial, and frequently very difficult, step encountered 
when using this technique involves identifying distinct features on the film (1) which 
appear on a sufficient number of frames and (2) the geodetic coordinates of which can be 
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determined. For many applications the user has little control over what Earth-related 
features will be photographed during flight and must instead work with what the film data 
offer. Many hours can be spent matching a series of features appearing on the film with 
the corresponding features on a topographical map so that their coordinates can be deter- 
mined. A feature whose coordinates have been identified is considered a landmark. 

Often, because of the flight trajectory and vehicle motions, a landmark may appear 
only on a relatively small number of frames. On the other hand, features appearing on 
the film which cannot be identified on a map may be visible for a relatively long period of 
time. Such nonpermanent features as clouds, ground vehicles, and ships might fit this 
category if they were slowly moving with respect to the camera frame rate (that is, quasi- 
stationary). The program has been modified to allow the user the flexibility to process 
several "trial" coordinates (latitude, longitude, and altitude) for a single feature in one 
computer run; this procedure allows for a quick iterative solution of its actual location. 
The user can accomplish this determination of the location by comparing vehicle Euler 
angle results obtained by using only known landmarks with those obtained by using both 
known landmarks and features. Solving for the coordinates of the feature readily follows 
and permits the feature to be used as a known landmark in other frames of the film where 
actual landmarks are not distinguishable. 

In applying this option, the iteration procedure is initiated by first estimating the 
location of the feature and running a series of latitudes about that estimation. By com- 
paring the resultant Euler angles with the values from the known landmarks, the best lati- 
tude is obtained. A series of longitudes is then processed at that latitude; as a result, a 
best estimate of longitude is obtained. The best estimate of altitude is established in a 
similar manner. With these new altitude and longitude values, the process can be 
repeated to define the latitude more accurately and, subsequently, the longitude and alti- 
tude parameters. The landmarks defined in this manner may be used to find the coordi- 
nates of other quasi-static features during other parts of the data period. Application of 
this option permits obtaining vehicle orientation data for time periods when no true land- 
marks are in view of the camera. 

As previously stated, a minimum of two landmarks is required for each frame. 
Theoretically, no limits are required on the maximum number but satisfactory results 
have been obtained by using from two to six landmarks, three landmarks being the pre- 
ferred number. The current computer program in use at Langley Research Center is 
set up with a six landmark maximum. 
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Basic data reduction ,- The basic data consist of film (image space) coordinates of 
the selected landmarks (or features to be converted to landmarks) with respect to the 
center of the frame. The basic data are punched on digital computer cards. 

By using the Gerber film reading system employed at the Langley Research Center, 
the basic data are obtained as follows: First, the processed film from the onboard 
vehicle camera is inserted into a film reading system capable of providing coordinates 
for either 16- mm, 35-mm, or 70-mm film. The film reader magnifies a 16-mm frame 
image to dimensions of approximately 26.7 cm by 37.3 cm (10.5 in. by 14.7 in.) and pro- 
vides coordinate readings to the nearest 0.00254 cm (0.001 in.) on the enlarged image. 
The film reading system is coupled with a visual- display electronic digitizer, a card 
punch, and an electronically operated typer (data typewriter). The card punch produces 
the basic data deck and the typer provides a record which can be used to determine 
rapidly whether any gross reading or film reading system errors exist. Then a manually 
controlled pair of crosshairs located on the film reader are consecutively alined on each 
corner (or other fiducial point) of the frame. These crosshairs also drive the digitizer. 
The coordinates and other reference information displayed on the digitizer are next 
recorded by both the card punch and typer by the user pressing a foot switch. For 
accuracy purposes, each reading is taken three times. The software is employed to 
average the three readings for each of the four corners of the frame and subsequently to 
determine the coordinates of the frame center. In the same manner the film coordinates 
of each landmark are obtained. 

PROGRAM DESCRIPTION 

The computer program (Program VOPD; Library No. A4424) used at the Langley 
Research Center is written in Fortran IV and is set up for a CDC 6600 series computer. 
For applications utilizing an average of three to four landmarks per frame, computer 
running time has averaged 0.17 second/frame. The basic data deck contains the film 
coordinates of the landmarks obtained as previously discussed. The control data 
include geodetic coordinates of the landmarks, vehicle trajectory data, and camera 
characteristics. 

Program VOPD is the main program. It calls subroutines TGTCAM and ALPBET 
for auxiliary calculations. A flow diagram for VOPD which indicates the basic calcula- 
tion procedure follows and a computer listing of the program is given in appendix B. 
Computer library routines MATINV and FTLUP are called by VOPD and are described 
in appendices C and D, respectively. 
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The following section gives a summary of the basic program operations and is 
followed by a more in-depth input-output discussion of the control data format and some 
sample program applications. 

General Operation 

The program operates under the following general steps: 

(1) The control data and basic data are read 

(2) The frame center and each landmark location on the frame is computed and 

stored on a disk for each frame of the basic data deck 
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(3) All landmarks are transformed from the geodetic to the geocentric coordinate 

system and stored in core 

(4) Each frame of data is then processed as follows: 

(a) The first frame is read from the disk and compared with the control data 

to determine whether it is to be ignored or whether any landmarks are 
to be ignored 

(b) The object space landmarks are converted to their image space coordi- 

nates and lens distortion corrections are computed 

(c) The camera Euler angles are solved and compared with the initial esti- 

mates from the control deck 

(d) Through an iterative process using the new Euler angles as a starting 

point, the final Euler angles are obtained when a sufficiently small 
change is detected 

(e) These camera angles are then converted to the vehicle-oriented values, 

and by using the meteorological and trajectory data of the control deck, 
the vehicle angle -of -at tack components are computed and printed 

(5) At the completion of all frames of the basic data deck, the program returns to 

the input group specified by the control deck for further processing. 

An initial estimate of camera Euler angles is necessary to start the iteration pro- 
cedure. For the first frame, the initial estimates should be in or adjacent to the quad- 
rant of the final solution to obtain satisfactory convergence within the 22 iterations pro- 
vided for in the program. For subsequent frames, initial estimates are taken as the 
results of the previous frame. 


Program Usage 

Basic data card format .- As previously discussed, the basic data consist of film 
coordinates of the selected landmarks with respect to the center of the frame. The basic 
data card should include, as a minimum, these data and the information specified in the 
program control-data deck. For program VOPD, this information includes landmark 
number and frame number. A separate card is required for each of the three readings 
of each point, whether it is a fiducial point or landmark. Four separate fiducial points 
are required for each frame. Thus, a frame with six landmarks would require 30 basic 
data computer cards. 

The basic data card should also contain the information required for identification 
and review, and also, as an aid for determination of possible reading or system errors. 
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The format used in the most recent Langley applications is listed in appendix E. Also 
included in appendix E is a typical sequence of steps employed in reading a frame. 

Control data format .- For the convenience of the potential user, program VOPD 
has been divided into distinctive input and output groups of control data cards. Detailed 
descriptions in computer terminology of the input for each group are included in appen- 
dix F. These descriptions include such information as column location, parameter sym- 
bols, parameter description, the number of cards in each group, parameter units, and 
input format. The following discussion, summarizing the information in each group, is 
in the order in which each group appears in the deck, and should be reviewed in conjunc- 
tion with the more detailed information given in appendix F. 

Group 1 contains a table of camera lens distortion compared with radial distance 
from the lens center. These values are used for, correcting the landmark image space 
coordinates and are obtained from camera calibration data. 

Group 2 is a title card to annotate each page of output. A "1" must be punched in 
column 1 to operate the carriage control on the printer and any descriptive information 
can be punched in columns 2 to 80. 

Group 3 consists primarily of data describing the conversion of camera Euler 
angles to vehicle Euler angles and is strictly a function of the relative orientation between 
the camera and vehicle axes. To solve for vehicle attitude angles accurately, it is essen- 
tial to know the exact orientation between the camera and vehicle axes. The numbers of 
values in the vehicle trajectory table (group 7) and wind table (group 4) are also required 
in group 3. The conversion data for camera to vehicle Euler angles are the coefficients 
of the equations defined in the input description (appendix F). Camera and vehicle Euler 
angles are shown in figures 1 and 2, respectively. 

Group 4 is a group of cards each containing an altitude and the corresponding hori- 
zontal wind components. The span of altitudes in this group must at least include the 
altitude span of the vehicle trajectory group (group 7). These data are required only for 
angle -of -attack and velocity calculations. 

Group 5 lists the number of landmarks in the basic data deck, the number of frames 
of data in the basic data deck, and the number of values in the time-frame table (group 6). 
The remaining spaces on this first card and all the spaces on the following cards in this 
group, as many as necessary, comprise a table consisting of the number of landmarks 
for each frame. These numbers are listed in the order corresponding to the frame 
sequence in the basic data deck. 

Group 6 consists of a set of cards each listing the frame number and an associated 
time. This group must at least include the frames in the basic data deck. For a constant 
frame rate camera, only two values are needed for this group. 
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Group 7 lists the vehicle trajectory data. The origin for this trajectory infor- 
mation is described in group 8. Each card of group 7 contains vehicle flight time, 
associated vehicle position (X, Y, and Z locations) and velocity, and altitude above 
mean sea level (MSL). This set of data must encompass the time span of the basic data 
deck. 

Group 8 contains the following data on one card: camera focal length; the origin of 
the vehicle trajectory coordinate system consisting of geodetic latitude, longitude, and 
the Earth radius to this origin; and the initial estimates of the camera Euler angles. 
These initial estimates are used to initiate the iteration on the Euler angles for the first 
frame of each basic data deck. Subsequent frames employ results from the previous 
frame to initiate Euler angle iterations. 

Group 9 lists the landmarks and their locations. This group includes all the land- 
marks read from the film using the same landmark identification numbers as used in the 
basic data deck. Geodetic latitude, longitude, and altitude for each landmark are listed 
on each card along with the landmark identification number. 

Group 10 indicates the number of frames from which specific landmarks are to be 
ignored. This value is, in effect, the number of group 11 cards. Landmarks are ignored 
when it is suspected that either their frame coordinates or Earth coordinates have been 
incorrectly identified. 

Group 11 consists of a number of cards as specified in group 10 which lists the 
frames and landmarks which are to be ignored for these frames. To ignore the entire 
frame, an option is also included as described in appendix E. 

Group 12 defines the following two parameters: First, the input group to which the 
program returns after processing all frames of the basic data deck. This information 
allows the user to stack jobs by returning to any input group except 11 and by adding 
behind the basic data deck only that input group and those following the one specified. 
Additional data in group 12 consist of the landmark number of any landmarks which are 
to be ignored in calculations for all frames of the basic data deck and eliminates the 
necessity for a long group 11 table. 

After these groups of cards, the basic data deck is inserted and appears only once 
in the deck setup. For subsequent processing, the basic data deck is automatically 
stored on the disk and need not be repeated. The information read from the basic data 
deck are the landmark number, its coordinates on the frame (x, y), and the frame num- 
ber. The landmark number is zero for the four fiducial (frame corners) points. 

Program output is listed by groups in appendix G and contains such information as 
parameter description, symbols, and units. 
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EXAMPLE APPLICATIONS 


For a better understanding of program application, several example cases taken 
from BLDT applications are presented with the corresponding program input and output 
format (for one frame only) presented in appendix H. Example case inputs are divided 
into the previously discussed input groups (appendix F) and may be referred to for clari- 
fication. For all cases, input group 12 provides instructions to the program to be imple- 
mented upon completion of the existing case. Program output consisting of a single 
frame for all cases is included collectively after the program input for all cases. 

Case 1.- This example represents an application where 24 frames are processed 
(only one frame is presented) and where no question exists concerning landmark loca- 
tion. All 12 input groups, as defined in appendix F, and the basic data deck are included. 
Three landmarks, identified as landmarks 65, 66, and 67, are used and their coordinates 
are listed in input group 9. As indicated in group 10, no landmarks are to be ignored. 
Group 12 provides instructions to the program upon completion of case 1. 

Case 2 .- Case 2 is a repeat of case 1 except that one landmark (landmark 65) was 
ignored because there was some question about its exact location on the film. The inputs 
for all groups, except group 12, are the same as those for case 1. 

By comparing the outputs of cases 1 and 2, a difference of about 1° can be seen 
for xfs and <p with little change in 9. By ignoring landmark 65, the residuals are 
slightly less for case 2 than for case 1. This result indicates that the results from 
case 2 may be more accurate than the results from case 1 and that the coordinates of 
landmark 65 may be in error. 

Case 3.- Case 3 is similar to case 2 except that the latitude of landmark 66 was 
changed from 32°50' to 32°50.3\ (See group 9.) The inputs for all groups, except 
groups 9, 10, and 12, are the same as those for cases 2 and 1. Case 3 illustrates a 
typical iteration to define the object space (Earth) coordinates of a previously unknown 
landmark to convert a feature on the film to a usable landmark. The results of this 
iteration show a slight improvement in residuals over case 2 and about a 1° change in 
and 9 . 

Case 4.- This case illustrates the combined approach of completely ignoring 
frame 55 and ignoring landmark 69 for several frames (frames 60 and 65) as instructed 
through input group 11. This case is independent of the previous three and requires a 
different set of input data. For this case, the output for frames 60 and 65 are presented. 
A computer listing of the general camera orientation method program (VOPD) is shown 
in appendix H. 
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CONCLUDING REMARKS 


Details of application of a previously devised photogrammetric method to deter- 
mine a time history of vehicle flight attitudes have been included in this paper. Emphasis 
has been placed on the techniques involved in reducing the raw photographic data to com- 
puter inputs in Fortran IV language and on the computer techniques and programs 
involved in obtaining vehicle flight attitude results. Also discussed are the major pro- 
gram modifications which allow faster data reduction and permit the user to determine 
the Earth-related coordinates of unknown or nonpermanent features appearing on the 
film. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., May 3, 1974. 
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APPENDIX A 


PHOTOGRAMMETRIC DETERMINATION OF CAMERA ORIENTATION 


In order to describe adequately how the camera orientation can be determined, it is 
necessary to define the coordinate systems in which observations are made. The coordi- 
nate system in which the position of observed points are known and in which the camera 
is oriented is called the object space. The coordinate system composed of the camera 
focal plane and focal axis and in which image coordinates are measured is the image 
space. (See ref. 3.) The relationship between a point in the object space and image 
space is 


X-Xp = C 

y - y p = c 


XjX + Mi Y + u 1 Z) 
AgX + Z j 

A^X + \x 2^ Y + 1^2 Z\ 
AgX + Z j 






(Al) 


The relationship in equations (Al) can be derived with the aid of figure 3. Let Q(X,Y,Z) 
be a point in the object space, then the image of Q will be l(Xj,Yj,Zj) in the image 
space. Since the origins of the two coordinate systems coincide for all practical pur- 
poses, the following transformation describes the coordinates of Q relative to the image 
space 


v 


A 1 /-t l ^1 


Y 

y i 

= K 

X 2^2 V 2 


Y 

_ z i. 


_ X 3^3^3_ 


Z 


where K is a constant of contraction. In equations (A2) dividing the first and second 
equations by the third equation removes the constant K and gives 


Xj A. jX + |i jY + v ^ Z 
Zj AgX + /igY + 

> 

Yj AgX + /igY 4- 

Zj AgX 4- M 3 Y 4- Z/gZ 


(A3) 
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APPENDIX A - Continued 



Figure J>.- Axis systems relating image space and object space. 

Letting 



and multiplying both sides of equations (A3) by c gives 


/XjX + \i ]Y + v^Z) 

X " Xp = VV c+ 'ti Y+ ''3 z / 

p ~ \ X 3 X + M 3 Y + V 3 Z , 


(A4) 


Equations (A4) are identical to equations (Al). 

For each known point in the object space there exist two equations (Al) relating the 
object space to the image space. The directional cosines of the 
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APPENDIX A — Continued 


image space axes relative to the object space axes can be expressed in terms of three 
angles a^, o^, 311(1 03. Consequently, these angles can be used to describe the orien- 
tation of the camera coordinate system relative to the object space coordinate system. 
To obtain the directional cosines, consider three successive rotations through the 
angles o^, and °3* By using figure 3 and imposing a constant of contraction for 
the camera, the following transformation is obtained: 


y - y r . 


= K 


LUO - Dili Uj 

cos a i sin cr^ - sin sin- cos a 3 
sin Dj cos $2 


sin cos - cos a ^ sin Og sin <7g 

sin sin - cos cr^ sin ag cos o g 

COS Cj cos O2 


cos 02 sin 


X 




X 

cos Og cos 


Y 

= K 

x 2 ll 2 v 2 


Y 

sin erg 


Z 


X3M3V3 


Z 


This equation is equivalent to equation (A2). 


The Least- Squares Solution 

The solution of equations (Al) contains six parameters which under theoretical 
conditions are constants. These parameters are Xp, y , c, Oj, o^, and Og of which 
5L, y , and c are measured independently for this experiment and are not unknowns in 
the solution. A solution of equations (Al) for cr^, 311(1 a 3 can be found with two 
properly chosen observations. Associated with x and y in equations (Al) are errors 
e and e. Since these errors exist, a computational method is needed which yields the 
best possible results with all the information available. The method of least squares 
which is described subsequently uses a minimum error criterion and has been used in 
the data reduction for this investigation. 

In general, equations (Al) with the associated errors can be written as 

*i = F (°j) + e i (j = 1, 2, 3) 

y i = F(oj)+ 6i (i = 1, . . ., n) 



where. F and F are nonlinear functions of Oj and in order to find a solution they 

must be linearized. Expanding equations (A5) in a Taylor’ s series about a nominal set 

a.: _ and dropping the higher order terms results in the following linear approximations: 
J? ^ 


A *i = *i - *i ,0 = b il( a l - a l,o) + b i2( CT 2 - 
A ^i = h - yi, 0 = b ii( a i - a i, 0 ) + b i2( a 2 - 


a 2,o) + b i3 (°3 - °3,o) + e i 
a 2,o) +b i3( a 3- a 3,o) + F i 


(A6) 
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APPENDIX A - Continued 


where 


b = !Ii 

il 9 (Tj 


9F. 


b. 0 = 


i2 - ac 2 
h 8F ‘ 

n i3 - 


°l,o a 2,o a 3,o 


a l,o ff 2,o cr 3,o 


a -4 (Tn a 0 

lj o 2,o 3 , o 


h aFi 

b il" 9 


h -!!i 
12 aa 2 


r 

13 8ct 3 


a l,o a 2,o a 3,o 


a l,o a 2,o CT 3,o 




Letting 

AC, i “ {°j - °j,o) 

equations (A6) can be put in the following form: 


- Act- + e - 
i] ] i 


aX i - 1 b i S 

j=l 

3 

A h " 1 b ij Ao j + h 
i=l 


(i =!,■.. n) 


(i = 1 n) 


> (A7) 


For further considerations the linear equations (A7) corresponding to the ith observation 
are expressed in matrix notation 


where 


y i = 






“l 


Act^ 

Aa 2 

_ 4tr 3j 


r 

'*4 

_ A \ 

B i = 

" b il 

b il 

b i2 

b i2 

b i3l 

b i3_ 

Act = 

•«-! 

e i 

A 


Vj = Bj Act + e A 


(A8) 
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APPENDIX A - Continued 


Then for n observations there are n matrix equations of the form of equation (A8) 
which may be written 

V = B Aa + e (A9) 

where 


V 1 


" B l" 


e l 

v 2 


B 2 


e 2 

• 

B = 

• 

e = 

. 

_ v n 


B n 


e n 


The problem may be restated: given V and B find the best estimate Aa 

for Aa. 

The best estimate Aa is the value of A a which minimizes the sum of the 
squares of the residuals e T e where 

e T e = (v - B Aa) T (v - B Aa) (A10) 


In order to minimize equation (A10), the first variation 6 with respect to Aa must 
vanish; that is, 


6(e T e) = e[(y - B Aa) T (v - B Aoj) = 0 
s(e T e) = -2 (v T - Aa T B T )B6 Aa = 0 J 


(All) 


Since 6 A a * 0, equations (All) can be satisfied if 
(v T - Aa T B T )B = 0 


or 

B T B Aa = B T V 

Solving for the estimate of Aa in equation (A12) gives 
Aa = (b T b)" 1 B T V 


(A12) 


(A13) 
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APPENDIX A - Continued 


A second necessary condition for equation (A10) to be a minimum is that the second 
variation with respect to Act be positive definite. Upon examination, the second varia- 
tion is 

6 2 (e T e) = 26 Act T B T B6 Act 

which is positive definite. Therefore, equation (A 13) is a valid expression for Act. 

Since equation (A13) is based on a linear approximation with nominal CTj q , Act 
can be used to find the best estimates <h. With the relationship a = a + Act, the value 

XV J U XV 

of A a which minimized equation (A10) leads to a new nominal m m n + Aa. This 
process implies an iterative procedure which continues until Act — 0 and the value 
of ctj 0 that leads to this result is the best estimate of ctj for crj. 

Partial Derivatives of Projection Equations 
For equations (A5) 

*i = F ( a j) + e i 0 = !* 2 > 3) 

Yi = F(ct.) + e i (1=1,..., n) 

The partial derivatives of F and F with respect to <jp ct 2 , and Og are as follows. 
Let 

p = A jX ,+ jY + 
q = XgX + pgY v 2 ^ 

r = AgX + PgY + 

then 
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APPENDIX A - Continued 



3F 




90^ 

r2 



cfJSLr- 

JL q \ 

9F 


d °2 ) 

9a 2 ' 

r 2 



cfJSLr- 


9F 

V°3 

9ct 3 / 

9a 3 

r 2 



JP = 

9(7^ 

JSU 

90-j 



JE = 
9<t 2 


_9q 

9a 2 


9r _ 
9o 2 

9a 3 


9q _ 
9a 3 


MjX - 

X 1 Y 






M 2 X - 

A 2 y 






M 3 x- 

x 3 y 






-^(x 

sin o j 

+ Y 

cos 

a l) 

- (sin a 2 sin 

"3) Z 

-"2( X 

sin 

+ Y 

cos 

a l) 

- (sin a 2 cos 


-u 3 (X 

sin aj 

+ Y 

cos 

a l) 

i + Z cos a 2 


X 2 X + 

4 2 y + 

"2 Z 





X 

1 

+ MjY 

+ i/ l 

z ) 





9r 

9a 3 


0 
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APPENDIX A - Concluded 


Coordinate Transformation 

Object space points and the position of the camera are initially identified in terms 
of geodetic latitude, longitude, and altitude above sea level. These data are obtained 
from maps of the photographed area and radar observations of the vehicle trajectory. In 
order to reference the data relative to the camera as described previously, the Earth- 
centered geocentric coordinates of both the object space points and the vehicle position 
are computed by 

X i = ( R i + h i) cos cos 
Y i = ( R i + h i) cos sin A-i 
\ = ( R i + h i) sin 

and 

^ - ll'35.6635" sin 2^ + 1.1731” sin 4^ - 0.0025” sin 6^ 

= 6378.388(0.998320047 + 0.001683494 cos 2^ - 0.000003549 cos 4$) 

where 

i = 1 for the vehicle position 

i = 2 . . . n + 1 for the object space points 

In the preceding discussion the X axis is in the equatorial plane pointing toward 
the Greenwich meridian, the Y axis 90° east in the equatorial plane, and the Z axis 
toward the north pole. A slight error is introduced by adding altitude to the Earth’s 
radius vector but this error is negligible for the accuracy desired of this system. 

The following transformation maps the geocentric coordinates into object space 
coordinates relative to the camera where the X axis points east, the Y axis points 
north, and the Z axis points toward the zenith. 


" x i + l~ 


-sin cos Aj 0 


_X i + l - X l _ 

Y i + 1 

= 

-cos A^ sin -sin sin cos 


Y i + 1 - Y 1 

>i_ 


cos Aj cos <f>j sin Aj cos sin 


5+1 ■ \ 


This transformation is obtained by a positive rotation about the Z axis through an 
angle of 90° + Aj followed by a position rotation about the new X axis through an 
angle of 90° - 
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APPENDIX B 


PROGRAM LISTING 

The program listing follows: 

PR0GR4M VOPO ( INPUT»OUTPUT»TAPE5=INPUT,TAPE6=OUTPUT»TAPE10,TaPE9) A 1 
DIMENSION WRR(2S), ISAVE(-JO), XSAVE(30), YSAVEI30), XRR<25). RR(20 A 2 
1) * DSR ( 20 > , OX Y ( 20 ) A 3 

DIMENSION XLAM ( 20 ) * XMU(?0). XNU(20>* R<2,6)» BS(6,2), XN(b,6), EP A A 
l S ( 2 , 2 0 > * C ( 6 , 1 ) » XMM ( 20 ) . XNNI20). QQ(20), nG( 20)» X(20). Y(20), S A 5 
2x ( 20 ) « 5 Y (20) ♦ IOENT (o) » CM<20>, CN(2u), IPIVOT(b), InDEX<6,2) A 6 

DIMENSION Ph 1 60 (10 01 , PH I UM ( 1 0 0 ) , PHlBS(lOO) ♦ XLOBD ( 10 0), XL08M<10 A 7 
10) A 8 

DIMENSION XL06S (100) * TX YZ T ( 20 , a ) * TGT T Ab ( 20 * 3 > * XYZTAB<3>» XX(20) A 9 
1* rr<20)» XNTEMP(b,6) , ANUNIT(bib) A 10 

DIMENSION TTH J ( 90 ) , TThJmI (mO) * TPH J ( 90 ) * TPH JMI ( 90 ) * TRJI90) A 11 

DIMENSION GEOUA (6,3) » GE0Ud(o,3) A 12 

DIMENSION NPTI200), IKRM(lOtl), TIMX(IOO). XXX(50)» YYYI50), ZZZ<50 A 13 
1) . XFkMUOO) A 19 

DIMENSION lournoo)* IFM(IOO)* JOUT(b.lOO) A lb 

COMMON /Ad/ Ill»lwV,TX(bO) » VVX (50 ) ,VVY (50) » VVZ (50 ) »H (50 ) , ALT 150) » V A lb 
: 1X0(50) * V YO < SO ) A 17 

COMMON /OIPCOS/ XLAM, AMU, XNU A IB 

DATA <kAO = S7.*;9577951 ) , (Pl=3. 1A1S926S) , (ETPKM=3280 .8336) ♦ (FACTOR=. A 19 
1713616A0020S) , (EkAUKM = b37B. 165) , <V = 0.0> A 20 

OATA (SFOC4L=.0015) , <SLAT=.083333J3333333) , (SXYZ=200. ) , (SRJ=200.) , A 21 
1 (SxPYP=.01) , (SAY=0. 1) , < VFTP<M=.00030M8006) , (RN=0.) A 22 

TwOPI=2.*PI A 23 

I TG T0= l A 2A 

OELF=0. A 2b 

1 READ (5.76) <XRR(I> ,RRR(I) ,1=1,25) A 26 

2 READ (5.77) A 28 

IF (EOF »5) 13,200 A 29 

200 j U m P = 0 A 31 

3 READ (5,78) 1 1 I , I WV , S 1 M, SI A , S2m » S2 A , S3M , S3 A A 32 

A READ (5,79) ( ALT ( 1) , VXU ( I ) , VYO ( I ) , 1=1 . 1 WV) A 3A 

5 READ (5,80) NT.NFR.nTFK, (NPT(L) ,L= 1.NFR) A 3b 

b READ (6,81) ( 1FRM ( 1 ) ,T1MX ( I ) , 1 = 1 ,NTFr() A 38 

DO 7 0 0 I = 1 , N TF k A AO 

TOO XF«M(t)=IFRM(l) A Al 

7 READ ( 5 » d'* ) ( T X ( I ) ,XXX( I ) , YYYI | ) ,ZZZ ( I ) ,VVX(1) ,VVY( I ) ,VVZ ( I ) .HI I ) , A a2 

11=1,111) A a 3 

3 READ (5,82) CC.PHI .THZ.KZ, ALPHAI .XOMEGI .XKAPI A A6 

PhITEm=PHI A A9 

T*CPhI=(2.*PHl)/R40 A 50 

FORPhI=2.»TWOPHI A 51 

SlxPnI=3.*TWOPHI A 52 

SInTP0=5IN(TW0PHI) A 53 

SINFP0=5In(F0PPHI) A 5A 

SINSPD=5IM(SIaPHI) A 55 

Ph I =Ph1-.1 9323986111 1 1°SINTP0* . 000325861 I'll 1 1 1 1*SINFPD»7 . 222222222 A 56 
1222E~07°SINS P 0 A 57 

PH I Z = P.H I A 58 

DPR=S7. 29577951 A 59 

PH IZ=°hIZ/QPR A 60 

THZ=THZ/OPR A 61 
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APPENDIX B - Continued 


9 D01 1 0 0 I JK=l *NT 

E£A0 ( 5 » a 3 ) itgt»thj.thjmi*phi j.phijmi ,rj 

RJ=kJo 1000./FTPKM 
7 TbJ 1 1 IGT)=Trtj 
TIfiJMl ( l TOT ) =TH JM I 
TPm J ( l IGT ) =PH I J 
TPfioM[ ( l T GT ) =PhI JMI 
1100 TkJ ( I TG t) =R J 

10 K£AO ( S * a 0 > NOUT 

IF (NOur.LT.l) GO TO 12 

11 F£AO <5»84) (IFM(l) * <JOUT(J»I> »J=1»6) ♦l=l»NOUT> 

DO 1 3'j I = 1 « NOUT 

DO 130J=1«6 

IF (JOUT(J»I).LT.l) JOUT 1 J* I ) =1 01 
130 COnTImj£ 

12 R£A0 (S*80) NGO* ( IOUT(I) , 1=1 ,NT) 

DO IS LL=l,NT 

IF ( tOuTlLU .LT.ll IOUT(lU=101 

15 CONTINUE 

IF < jUmP.GT.0) GO TO 20 
DO H N= 1 * NF P 
IUP = NP r (N) <*3* 1 2 

Ft ao (s«as) (i save m »xsave<i> *ysaveii> «ifram,i=i»iup) 

vpc=o. 

xpc=o. 

00 16 1 = 1 » 12 
xpc=xpc*xsave ( n / 12 . 

16 YPC=rPC*rSAvE(i)/i2. 

I UP = NP TIN) 

11 = 12 

DO 18 I=l*IUP 

>SAv£m=a. 

YSAve<n=o. i 

00 1/ J= 1,3 l 

h = ii>i | 

XSAVE ( I * =XSAVE ( II *XSAVE( II ) /3. : 

17 YSA VE < n’ = rSAVE ( I) *YSAVE ( 1 1) /3. ' 

18 1 SA vE ( 1 1 = I SA VE ( 1 1 ) 

JUMP=1 

20 k£x I NO 0 
REA0FG=0 
IF OOP G = 1 
KK = 0 

aupmap=alphAI 
xOmEGP=XO m EGI 
X K App = aK AP I 

21 KK=KK*t 

IF (KfC.LE.NFR) GO TO 22 

ENO FILE 10 
* rt IT t. (b»77) 

GO TO ( 1 1 2 » 3 * 4 >5 » 6 * 7 » 8 » 9 * 1 0 « 1 1 » 1 2 » 1 3 ) » NGO 

22 NJ=NPT<<K) 

RE AO (v f d7) IFPAM,XPC.YPC 
XFRA*f.= (FRAM 

CALL F1LUP ( XFRAM, TI * 1 *NTFH« XFHMtTIMX) 

CALL FILUP (TIfXYZ,l*III,TX.AXX) 

XYZTA-31 1)=XYZ 

CALL FTLUP (TI *AYZ»1 *11 I.TX* YYY) 

XY2TA8 (2) =XY2 

CALL PTLUP <TI*XYZ*1»IU»TX,ZZZ) 

AYZ TAJ ( )) = X Y Z 


A 62 
A 63 

A 65 
A 66 
A 67 
A 68 

A 71 
A 72 
A 73 
A 7* 
A 75 
A 76 
A 77 
A 79 
A 80 
A bl 
A 82 
A 84 
A 65 
A 88 
A 87 
A 88 
A 89 
A 90 
A 91 
A 92 
A 93 
A 94 
A 95 
A 96 
A 97 
A 98 
A 99 
A 100 
A 101 
A 102 
A 104 
A 105 
A 106 
A 107 
A 108 
A 109 
A 110 
A 111 
A 112 
A 113 
A 114 
A lib 
A 117 
A 118 
A 119 
A 120 
A 121 
A 122 
A 123 
A 124 
A 125 
A 126 
A 127 
A 128 
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APPENDIX B - Continued 


IF (nOuT.LT.I) GO TO 25 A 129 

DO 29 I = l • NOUT A 130 

IF ( IFkAm.n-I. IFM(I) ) GO TO 24 A 131 

IF (JOUTU.I) .NE.100) GO TO 25 A 132 

DO 23 L=1»NJ A 133 

23 RE AO (9.67) ITOT.XCO.YCO A 139 

• RITE (6 . 88 ) IFRAM A 1313 

GO TO 21 A 136 

29 CONTINUE A 137 

25 INJ=0 A 1.3a 

NK=0 A 139 

WRITE (6.77) A 190 

WRITE (o*89) A 1 h1 

DO 31 JJ=l»NJ A 1^2 

PE AO (9.87) ITGT.XCO.YCO A 193 

IF (N0UT.LT.1) GO TO 26 A 199 

DO 27 1 = 1. NOUT A i<*s 

IF (IFkAM.NE.IFM(I)) GO TO 27 A i A 6 

DO 26 m=1 .6 a 197 

IF (ITGT.inE.JOOT(M.I) ) GO TO 26 A 198 

NK=NK*1 A 199 

IF ( (NJ-NK) .LT.2) GO TO 28 A 150 

WRITE (o*9 0 ) I T G T A151 

GO TO 31 A 152 

26 CONTINUE A 153 

27 CONTINUE A 159 

28 00 29 I=1.NT A 155 

IF ( I OUT ( I > .NE.ITGT) GO TO 29 A 156 

NK=nK. I A 157 

IF ( (NJ-NK) .LT.2) GO TO 30 A 158 

WRITE (6.90) ITGT f A 159 

GO TO 31 A 160 

29 CONTINUE A lbl 

30 INJ=INJ*1 A lb2 

X ( l N J ) = XCO a 163 

Y ( I N J) = 7C0 A 

ThJ=TThJ( ITGT) A lb5 

ThJM(=TThJMI (I TGT) A 166 

PHI J=TPHJ(ITGT> A i b 7 

Prlt JMl=TPnJMI ( (TGT) a 1 08 

RJ=IRJ([TGT) A169 

IDEnT (INJ)=.TTOT A 1 73 

Tri J= Tri J. TmJM ( /60 . A 179 

GEOO A ( I N J . 1 ) =TH J A 175 

TH J=Th J*Rn*SL A T a 1 7 (, 

GEOOd ( 1 NJ* 1 ) = ThJ A 177 

PHI J=PhI J*PH l JHI/60. a 178 

GEOOA ( 1 9J.2) =PHI J A 179 

PHlJ=PhIJ.KN*SLAT a 1 80 

GE003 ( INJ.2) =PHI J A lMl 

PHI J5A=RHI J ,. A 1h2 

GEOOA ( INJ.3) =wj A l8 3 

RJ=RJ*KN°SRJ A 189 

GEG06 ( INJ.3) =RJ a 1^5 

RJ=RJ/F TPKM.ERADKM A 

Ph I =P-i I J A 187 

T*0Phl=(2.«PHl )/RAO A 188 

FORPr-iI-rP.^TWOPHl A189 

SI XPr-!l=3.*T*0PHl A 190 

PHI =Prt I TE M A 191 

SInTP0 = 51N(T>/0PHI) a 192 

SINFP0=5IN(F0«PHI) A 193 

SINSP0 = SIN(SIX Ph I ) A 1 9 9 
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APPENDIX B — Continued 


31 

32 


33 


3** 

3b 

3b 

37 

38 

39 


FHJ J=s-il J-. 193239«ol 11 11 *blNTPO+. 00003258611 1 1 1 1 1 1 « S I NFPD* 7 . 222222 
1222222E-07«5INSPO 
1GT Tfl-j < INJ, 1 1 = TnJ 
TGTl4-3(I.NJ»2)=PhIJ 
TO T T 4 3 ( INJ.3) = RJ 
1 T G r P = U 
NJ=INJ 
N=,\J 

CCSAv£=0. 

I 809 = 0 

KUnFLG=1. 

IkUnC0=10 
hNS*Y2 = Ki\l«SXYZ 

at z ra-m >=xrzn«m ♦rnsxyz 

AY l I AS (2) = X Y Z f AB (2) ♦UNSXYZ 
xrZTAB(3)sXY2TAHr3> ♦ KfJSXYZ 
TAY2T C 1 * 1 ) =T I 

T X Y Z T ( 1 »21=AYZTA6(i)*VFTPXM 
TayZT < 1 »3> = XYZTA8 (2) *VF TP m 3 
TAYZT (1 ,9) =XYZTAo{3) ° VFTPKM 

CALL TGTCAM (0) »NJfRZ,PHlZ.THZ,TXYZT,TGTTA8> 

00 33 1=1. nJ 

TGI r A3 ( I ,3) =TGTTA8 ( I . 3 ) 1 PKM/1 GOO . 

tfAITE ( b , 9 I ) lOENTU) , (TGTTAd(IiJ) ,J=1,3) . (GE0DA<I,J) ,J=2,3> »XLAM( 
II) iO’OIII * XNIJ ( I ) 

IGT 1 A6 < l » 3) =TGTTAd (I ,3) /F f PKM<* 1 00 0 . 

CC0*IG=CC 
CC=CC*wM*3F0CAL 
ccl A S T =cc 

if (ccsa^e.gt.o.) cc=savf 
CC1=CC I 

ALPbA V = Al!.PHAP 

0M£ 3 1 v=X0MEGP 
JHASA \1 = AKAPP 
W*[1E (0.92) 

WP I I c (6.93) (XYZTABU) ,1 = 1.3) .XPC.YPC.CC 

APC=*PC*PM»SXPYP i 

YPC=YPC»PN*5XPYP 

APCsxPC^PACTOP 

YPC=YPC a FACTOR 

xPC3T0=XPC 

YPCbTO=YPC 

APCl=XPC 

YPC l = YPC 

APC 5 A v = APC 

YPC 5 A V = YPC 

TwCn=2.»n 

00 39 1=1, N 

AX ( 1 ) =X ( I) 

Y Y < 1 ) = r ( I ) 

DO 3b l = 1 . N 
A ( t » =X ( 1) ♦«V*SXY 
00 36 1=1, N 
Y (I ) =Y ( I ) ♦PN^bXY 
WRITE (6,95) 

00 3/ 1 = 1, N 


A 195 
A lvo 
A 19 7 
A 198 
A 199 

A 201 
A 202 
A 209 

A 207 
A 208 
A 209 
A 210 
A 211 
A 212 
A 213 
A 219 
A 215 
A 21b 
A 217 
A 218 

A 219 
A 220 

A 221 
A 222 
A 223 
A 229 
A 225 
A 22b 
A 22 7 
A 228 
A 229 
A 230 
A 231 
A 232 
A 233 
A 239 
A 235 
A 25b 
A 237 
A 238 
A 239 
A 290 
A 291 
A 292 
A 293 

A 299 
A 29b 
A 29b 
A 29 7 
A 298 
A 299 
A 2b0 


hfUTE (6,99) IOENT(I) »X(I) »Y(l) 

00 38 I = 1 , N 

Mll=AU)*FACTOR 

Y(I)=Y(1)*FACT09 

IF (IRUN) 92,92,39 

GO TO (9b, 90), IFOCFG 


A 251 
A 2b2 
A 2b3 
A 289 
A 2b5 
A 2b6 


27 



APPENDIX B - Continued 


4 C CC=CC*OELF A 2!?7. 

APC=XPC5AV A 256 

YPC= YPCSAV A 259 

APCl=xPCSAV A 2b0 

YPC1=YPCSAV A 261 

DO 41 I = 1 * N A 262 

A(I)=XX(l) A 2b3 

41 Y ( I ) =Y7 ( I ) A 264 

42 IHUiN=I4UN*l A 265 

IF ( IwUNCO-IRUN) 43,45,45 A 26b 

43 IF (4UNFLG) 44,44,32 A 267 

44 RUNFLG=1. A 268 

CC=CC1 A 2b9 

DELF=-0ELF A 270 

IRUN=0 A 2 71 

GO 10 40 A 272 

45 CONTINUE A 273 

I CA 6E = l A 274 

lCA 3 e=lCASE*l A 275 

ALPi-i A = ALPS AV A 27b 

AOm£G = 0“1ESAV A 277 

XKAP=XKASAV A 2 78 

ALP«A4=ALPHA/RA0 A 279 

OM£GAP=<Om EG/RAO A 260 

xkapr=xkap/rao A 26 1 

1L I NE = 0 A 262 

SlC-X=l.O A 263 

S I G 7= 1 • 3 A 264 

NOP 1 = 0 A 265 

NOP2-0 A 266 

£P0C = AO . 00 A 267 

iccu\r=<: a 266 

1 TE" = 1 A 269 

JJ=1 A 290 

SUM=0 • 0 A 291 

46 CONTINUE A 292 

IF (ICA5E-1) 48*46*47 A 293 

47 CONTINUE A 294 

IF CITtR-22) 50.48*48 A 295 

48 SIG1=5M*ALPHA*S1 A A 296 

5IG2=S2M*X0MEG*S2A A 297 

SIG3=S3'4*XKAP*33A A 298 

IF 1SIG1.LT. 0.) SIGl=SlGl*3oO. A 299 

IF (5IG1.GT.36U.) S I G l = 5 1 G 1 -360 . A 300 

IF (SIG2.LT. -130.) S(G2=Sl02Ob0. A 301 

IF (SIG2.GT.l60.) SIG2=SIG2-3bO. ' A 302 

IF (SIGi.LT.O.) SIG3=SIG3«360. A 303 

IF (SlG3.GT.3bO.) Sl31=SIG3-3bO. A 304 

WRITE (5*9 6 ) A 305 

DO 49 I = 1 * N A 306 

49 *RITE (6.97) Ii)ENT < I ) * X ( I ) » Y ( 1 ) *RR < I ) *0RW ( I ) A 30 7 

CALL A|_PdET ( T I » ALF A . MET A * E TA * C VP * S I G2 * S IG 1 . SI 63) A 308 

WRITE (6*96) A 309 

WRITE (6*99) IFRAM* TI * ALPHA « XOMEG * XKAP * XPC * YPC.SIG1 *SIG2*SIG3* ALFA A 310 

l.aETA.ETA A 311 

WRITE (10) TI,3IG1»SIG2*SIGJ. ALFA, bETA, ETA A 312 

50 CONTINUE A 313 

SALPmA = SIN(AI_PHAR) A 314 

CALPHA=C05(ALPHAR) A 315 

SOmEGA=SIN (OMEGAR) A 316 

COMEG A=C0S(0MEGAR) A 317 

SKAPPA=SIN(XKAPR) A 318 

C*APPA=COS(XKAPR> A 319 
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AAC=-CALPhA"CXAPPA-SALPHA*SOMEGA"SKAPPA A 320 

6aC=SALPrtA*CKADPA-CALPHA*S0MEGA*SKAPPA A 321 

ccc=comega»skappa a 322 

A ACP-C ALPrtA "S nAPPA-SALPHA & SOMEGa*CKAPPA A 323 

edCP = -SALPHAOiiKAPPA-CALPHA*S0M£GA*CKAPPA A 324 

CCCP = CCM£GA»C*A=>PA A 325 

D0C=SAL?nA»C0N;tGA A 32b 

EEC=C Ai_Pn A*C0MEGA A 327 

FEC = S0~(EGA A 32t) 

IF ( I CASE- 1 ) 52,S2,51 A 32V 

51 CONTINUE a 330 

IF (ITEN-22) 53*52*52 A 331 

52 CONTINUE a 332 

WHITE (0,100) A 333 

53 CONTINUE a 334 

j=0 A 335 

54 J=J*1 a 33b 

/MC = AAC'»XLAm( J) *bHC*AMU(J) ♦CCC°XNU(J) A 337 

XNCsAACP^XLAMf J) ♦driCP°XMU(J) ♦CCCP*XNU(J) A 33B 

CC = OOC°XLAN ( J) ♦ EEC*XMiJ< J) ♦FFC*XNU< J) A 339 

xm,v ( j) = AMC/CC a 340 

XNN( J» =XNC/GC . A 341 

CO(J)=CC/OC A 342 

AA1=5 )3T< U<J1-APCST0)""2*(Y(J>-YPCST0)"*2> A 343 

HK(Jt=XXl A 344 

CALL FfLUP (XaV,0ELH»1*25«XHH,HHH) A 345 

0£lh = -uELK/1000. A 340 

OKH ( J) = OcLK A 347 

OISfOP=l.-OELH/XXl A 34b 

EPSIl . J) =OISTOKMX< J>~APC)-CC«XMM< J) A 349 

fcPS<2,J>=L)lSrON<UY(J)-7PC)-CC<»XNN(U) A 350 

0*Y ( J1 = 50«T <EP5< 1 , J) "»2*FP5(2, J) "<*2> A 351 

CM ( J) =CC"AM4 ( J) A 352 

CM J) =CC»XNN ( J) A 353 

IF UCASE-1) 56,56*55 A 354 

55 CONTINUE a 355 

IF (I TEH-22) 57,56,56 A 356 

56 CONTINUE A 357 

WHITE (6,101) EPS ( 1 , J) ,tPo (2, J) » IOENT I J) A 35b 

57 IF (ITEH.LE.l) SLM=SUM+EPS< 1 , J) "EPS (1 » J) *EPS (2* J) "EPS (2* J) A 359 

IF (N-J) 5B,5b,54 A 3bO 

5b CONTINUE A 361 

1L INE= IL INE*6*N a 3o2 

IF (ILINE.GE.45) ILINE=0 A 3o3 

JJ=t A Jo4 

IF (ITEH.LE.l) XXN=N A 3o5 

UO 62 J=1,N A 3bb 

IF (J-nU(JJ)) 50,59, oO A Jo 7 

59 JJ=JJ»I A 3ob 

GO TO o2 A Jo9 

60 IF ( 4 3S (EPS ( 1 , Jl 1 -</) 61,61,03 A 370 

61 IF ( MS IEP6 ( 2 , J) ) - V ) 02,62,b3 A 371 

62 CONTINUE a 372 

WHITE (o* 102) ( (EPSTM, MM) »M=1 ,2) ,MP=1 ,N) A 373 

WHITE (a, 103) A 374 

GO TO 21 A 375 

63 GO 04 1=1,3 A 370 

C ( I * l ) =Ui 0 A 377 

00 04 N=l*3 a 37b 

XN ( I ,<) =0.0 A 379 

64 Con I l NUE A 3d0 

JJ=l A 3»1 

00 1C J = 1 » N A 3b2 

1 F ( J-nG ( J J) ) 66* b5»bb A3b3 
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oS JJ=JJ*l A 364 

GO 10 70 A 3b5 

t>6 FFl=aac-DOC*X*M(j) A 386 

FF2=86C-E£C<*xmm< J) A 3b7 

FF3=CCC-KFC*X.mM<J) A 386 

F31=aaCP-DOC*aN'\M J) A 389 

Fd2=63CP-EEC<»Ai\!N<J> A 390 

F63=CCCP-FFC*XNN< J) A 391 

GG1 :xm« ( J) *CO*FGA*FFC»SKAPPA A 392 

Gbl=XNN( J) «CO m£GA*FFC«CKAPPa A 393 

Hrtl = x(_A*i ( J) “SALPH A ♦ AHU ( i ) •'(J AlPHA A 39** 

S1=-FF2»aLAm(J) ♦FK1*XMU(J) A 39S 

S2=FF 3°riril *GG I “AnU ( J) A 396 

Sd l =-Fb2»AI_AH ( j) *Fbl*AMU< J» A 397 

S3?=F33 tt MHl *Gdl °AiJLM J) A 396 

b ( I . 1 ) =81«00< J) A 399 

H<l*2)=52‘»QJ<J> ' A 400 

b ( I » 3) = -CC*xni\i ( J) A 401 

H ( 2* It = 561 *00 1 J) ' A <*02 

h (etc) = '>b2<*00< J) A AO J 

B (2*3) =CC"*mm< J) A 404 

UO o7 I =1 t3 A 40b 

b’5 ( i 1 1 ) =b< 1 1 1 > A 406 

bb ( 1 * ?> =6 <2* I ) A 40 7 

67 COMIMOE A 406 

00 69 1=1*3 A 409 

CO 63 L= 1 * 2 A 410 

C<I.I)=C(I»1>+6S<1»L>*<-EPS(l*J>> A 411 

68 CONTINUE A 412 

00 69 K= 1 1 3 A 413 

DO 69 L= 1 * 2 A 414 

XN(1 .<<)=XN(I,K> ♦dS(ItL)*b(LtK) A 415 

XNTEMCM t»KJ=XN(I*K) ♦rfS{I*L)«d(L»K) A 416 

89 CONTINUE A 417 

70 CONTINUE A 416 

J1 =3 A 419 

- J2= 1 A 420 

CALL MAT I N V ( AN ( 1 » 1 > tJl.tCTTtl) *J2*U£TERM»IPIV0T* INDEX 1 6 1 1 SCALE) A 421 

00 71 1=1*3 A 422 

00 71 J= 1 * 3 A 423 

71 ANUNI T ( I * J) =0. A 424 

00 72 1 = 1 * 3 A 425 

DO 72 K=1 *3 A 42o 

00 72 J= 1 * 3 A 427 

72 XNUN[ T ( I * i< ) = XN ( ( tJ) oXNTEmMIJjK) ♦ XNUN IT ( I * K ) A 426 

ALPnAh! = ALFbAK + C( 1 t 1 ) A 429 

0m£G4R=0mEGAR*C<2* 1) A 430 

XMPR = AKAPR*C (3* l) A 431 

IF ( XXAPk.LT. 0. ) XKAPK = X*APr + T!«OPI A 432 

IF (ALPmAP.LT.O.) ALPHAH-=ALPHAR + TW0PI A 433 

IF (GMEGAR.LT.-Pl 1 OMtGAP=OM£GAP*T*NOPI A 434 

IF (xkaph.GT.TvOPI) AXAPR=XKAPR-TwCPI A 435 

IF < AtPbAP.(?T.T^OPI ) ALPhAR=ALPHAR-TwOPI . A 436 

IF I0M£GA4.OT.PI) OMEOAH=CMEvjAK-T/(OPI A 437 

Al_PrtA = ALPHApoRAU A 438 

X0MEG=0MEGA««4AU A 439 

XKap=XKAPP*PAO A 440 

IT£w=I TEK*1 A 441 

ICOUNT= ICOU.'JT + l A 442 

IF ( I COJN T-2 1 ) 46 *,46 *73 A 443 

73 SU*=0. ' A 444 

00 74 J= 1 * N A 445 
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74 

SO,- = 5 cM»E PS ( 1 * J) »»2»EPS (2, J> **2 


A 

446 


SUV=SQhT 15UM/TWUN) 


A 

A A 7 


wri rt 

( ij » 10a) mJM 


A 

448 


ALPH4P 

= ALPrlA 


A 

449 


XO M EGP 

= ROMEO 


A 

450 


XK4PP = 

XKAP 


A 

451 


IF (AUPHAP.LT.'I..0H.ALPHAP.GT.TWUP1) alphap=alphai 


A 

452 


IF ( XK4PP.LT.0. .0K.XKAHP.6T. TWOPI ) XKAPP=XKAPI 


• A 

453 


IF (xoi'EoP.lT. -PI.OP.XOMKGP.GT. pl) XUMEGP=XOMEGI 


A 

454 


GO TO 

(21 .AO ) t IFOCFG 


A 

455 

13 

WR I TE 

(rj.105) 


A 

45b 


STOP 



A 

457 

7b 

F OP*a T 

(2F10.3) 


A 

459 

77 

FOSoa T 

(60H 


A 

460 


1 



A 

46 1 

78 

FORMAT 

( 2 1 5 , 6r 5. 0 ) 


A 

4b2 

79 

FORMAT 

(3F10.1) 


A 

463 

a; 

F 0 A M A T 

(1615) 


A 

464 

81 

FORMA f 

(115.1F10.1) 


A 

465 

82 

F Ohm AT 

(bF 10.1) 


A 

466 

83 

FOfioA i 

(1 l 10.5F10.2) 


A 

4b 7 

84 

format 

(715) 


A 

4o8 

05 

FOP o A r 

(115/(715) ) 


A 

469 

80 

FOhoA j 

(8A» 1 12. 3GX,1F7.3,3X.1F7.3,6X,1I4) 


A 

470 

07 

FORMAT 

(1 l 5 . 2t 1 3. 7) 


A 

471 

08 

f op" a r 

(////6x, “FRAME*. 1 14, “IGNORED bY USERS DIRECTION*) 

A 

472 

89 

FOkmA T 

(/I X,*L-MAHK LONGlTU lE lat-gocn PAD £R KM 

LAT-GODT ALTir 

A 

473 


luot-p r 

XLAo XmU ANU*) 


A 

474 

90 

format 

(1 I5.2A. “IGNORED Hf UnERS DIRECTION*) 


A 

4 75 

91 

F Ohm AT 

(l 1 5 , 2 a » 2F 10. b, 1 F 1 (> . J « IF IU.5.1F1 1.J.3F10.7) 


A 

4 7b 

92 

fopoa r ( ino,»vEHfCLF a r z 

XPC YPC 

A 

477 


1 

CC“> 


A 

4 7 71 

93 

F OooA T 

( 7X.3F 10.2, JF10.5) 


A 

4 7b 

94 

forma r 

(If 4. 2F 15.6) 


A 

4 79 

95 

FOAM A y 

( 1h0.32hL-MAhK XCO FINAL YCO FINAL) 


A 

4b0 

9o 

FGft.oA 1 

( InO , *L-MAkN*.6X. l MX , 1 <, a • 1 rt Y » 1 1 X.6HH UIST.4X 

» 7H0ELTA R) 

A 

481 

97 

FORMA r 

(10 , I3.4F15.6) 


A 

482 

9h 

FORMAT (In0i» FWAMt TIME ->IGl SIG2 SIG3 

XP YP 

A 

463 


lPbt 

ToETa PHI ALPHA HtTA ETA*) 


A 

4831 

99 

FOPmA f 

( 115, Ah 1.3.2F7. J«?F IT. J, 1FH. 3, 3F8.3, 1F9.3) 


A 

464 

loo 

fo^ma r 

(lril», 47 H RESIDUALS RMI) RY(I) 

L-MARK) 

A 

485 

U1 

f orma r 

(2Xf2Flh.b,lllo) 


A 

486 

102 

F 0 9 o A T 

(//2E21 .5) 


A 

487 

103 

FORMAT 

( /// 1 l H ENO OF JOB) 


A 

488 

U4 

FORMAT 

(1H0,8X,*STANUARU DEVIATION OF RE S I DUALS* » 1 F 1 3 . 5) 

A 

489 

105 

FORMAT 

(//» 1 X,*ENJ OF case*) 


A 

490 


ENO 



A 

491- 


subroutine tgtcam (nc.nj.rz.phiz.thz.txyzt.tgttab) 

8 

1 

DIMENSION A ( 3 » 3) » XYZTAri(J»l), X5YSZS(3,1), TGTTAB (20.3) ♦ 8(3.3). 

8 

2 

ID A Y Z J ( 3, 1 ) * OxYZJP(3.l), AMV (3,1) , AMVB(3.1), OEL JPT (100.3) » TXYZT 

8 

3 

2(20.4) 

8 

4 

0 1 8c N 5 1 ON «TAb(20), AMUTA6(20) , VTA8(20) 

8 

5 

COmmOn /DIPCOS/ ATAH.Ar-iUTAb.vTAH 

B 

6 

OPR=57. 29577951 

B 

7 

C0SPZ=C0S(PHI2) 

B 

a 

CO$ThZ=COS(TH7) 

B 

9 

SInP7=S(N(Ph17) 

B 

10 

SINTH7=SIN(TH2) 

B 

11 

XZ=9Z*C08PZ»C05ThZ 

B 

12 
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YZ = -,Z»C05PZ*SINTHZ B 13 

Z2=4Z«5lNPZ B 14 

ASY5Z5 < 1 * 11 =XZ B 15 

ASy 5ZS <2,11 =YZ B lo 

XSY 5ZS <3» 1 ) -'LL B 17 

A { 1 . 1 » =-S(nThZ B IB 

All *21 =-Sl,\PZ‘ > COSTrtZ B 19 

A ( 1 *31 = »COB |J Z 0 COSTHZ B 20 

a<z,h=costhz b 2 i 

A (2*21 = -S I p»PZ*SINTHZ B 22 

A (2*31 =SInThZ b COSPZ B 2J 

A ( 3 ♦ l ) = 0 . B 24 

A(3,2l=CObPZ B 25 

A { 3*31 = :sl NPZ B 2t>- 

1NC=0 B 27 

1 1 NC = I NC ♦ l B 28 

Tl=r<YZT(lNC»t> B 29 

AYZ Ta-i ( 1 , l ) =IAYZT < INC, 2) B 30 

XYZfA-j(2«l)=TxYZT(INC*3) H31 

A YZ ras < 3. 1 ) =TaYZT ( INC»4) B 32 

TGTPLG=0. B 33 

INJ=0 B 34 

2 INJ=1NJ*1 B 35 

THJ0=TGTT4B<I,9J»1'> B 36 

Ph I JO= r 0 r T Ad ( I N J ♦ 2 ) B 37 

kJ=T3rrAB(INJ,3) B 38 

Thj=ThjJ/OPs B 39 

PhI J = Pri[ JO/OPk B 40 

IF ( TjfFLGl 3,3,5 B 4l 

3 COMInuF B 42 

IGTFl3=1. 8 43 

00 4 1*1 ,3 8 44 

DO 4 j=l ,3 B 45 

4 ASYSZb ( 1 , I ) =A ( I , J) “XYZTAb ( J*1)*XSYSZS< 1, 1) B 46 

XS=XS YSZ5 (1,1) B 47 

Y5=A5Y5ZS (2, 1) B 48 

ZS=x5Y5Zb(3,l) B 49 

T£MP=X5<»*2*YS«*2 B 50 

TEvP=S09r (TEMPI 8 51 

AL=a7aN2(ZS,TEMP) 8 52 

A LP r = AU "DPP 8 53 

T = ATa N z-(YS,x5) 8 54 

TPT = T «l)Pk 8 55 

TtMP*A5»<*2* YS‘‘ 0 2*ZS < ‘ 0 2 B 56 

P=50hT (TE^P) 8 57 

SIMT = SH(T1 8 5d 

COSI=CO->(T) B 59 

SIN8L*SIN ( AL> 8 60 

C05 A(_ = C'.)b < A L 1 fa 61 

fa ( l , 1) *-31 nT 8 62 

fall *21 = CUbT 8 o3 

fa ( 1 ,3) =ii. B 64 

fa (2, I) =-SINAl°CObT B 65 

fa(2.2l =-5INAL»blNr B 66 

b ( 2 , 3 ) -COSAL B 67 

8(3, I* =COsAl*COST 8 68 

bI3.2l=.ilNT»C05AL B 69 

fa ( 3, 3) =b INAL 8 70 

5 CONTINUE B 71 

COSPj=CU5(PnI j) 8 72 

C05Thj=CUS(ThjT 8 73 
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SIIm0J=SIN<Ph! J) 8 79 

SlNTnJ=SIN(TMj> 8 75 

XJ=kj*COSPJ»CUSTHJ 8 76 

YJ=RJ»COSPJ»SINTHJ 8 77 

2J=RJ*S1nPJ 8 78 

D£Laj=.VJ-XS 8 79 

OELYj=YJ-YS 8 '80 

G£LZJ=ZJ-ZS 8 81 

OELJPT ( 1NJ» 1 ) =OELXJ B 82 

D£LJ D T ( INJ»2) =0ELYJ 8 83 

0£ L J° T ( 1 N J » 3) =DELZ J 8 89 

T£mP = 3£LxJ*»2*OELYJ»*2*OELZJ**2 8 85 

CELRj=S )«T ( TEMP) 8 8b 

AKv ( l , l 1 = UELX J/OLLKJ 8 87 

AMy (2, l) = i)r.LY j/OELKJ 8 88 

AMV (3, l ) =O lLZJ/OELRJ 8 By 

00 o 1 = 1*3 8 90 

Amv6 ( I * 1 ) =0. 8 91 

GO 6 J= l * 3 8 92 

6 AMv8(1»1)=8(I»J)°AMV(J*1)*A-ivB(I*1) 8 93 

OXY Z J ( 1 , 1 ) = Jt L X J 8 99 

DxyZJ( 2,1>=0ELYJ 8 9b 

DxyZj(3»1)=0ELZJ 8 96 

DO 7 1=1*3 8 97 

DxyZjp(I,I)=0. B 98 

00 7 J= l , 3 8 99 

7 DA Y Z J J ( 1 * l ) =b ( I * J) tt UAYZJ(J» 1) ♦UXYZJQ1I * 1 ) 8 100 

ATA6(I<sJ)=AmV8(1»1) 8 101 

ArtUT A3 ( 1NJ» =A«.Vrt (2» 1 ) 8 102 

VTA8(InJ) = Am\/o(3*1) 8 103 

IF (NJ-INJ) 8.8*2 8 109 

8 CONTINUE 8 105 

DJPIFG=J. 8 106 

D JP T P G= 1 • 8 107 

IF (OjJfFO) 9,9,11 8 108 

9 hKI IE (o,l9) 8 109 

00 10 I N J - l ♦ N J B 110 

10 WRITE (0,13) (JELJPT < INJ,i\) ,K=1,J) B 111 

n continue h 112 

IF (NC-INC) lc.12,1 8 113 

12 RETURN 8 119 

C 8 115 

13 F OR T (lrt0*5x, 7E15.7) B 116 

19 FORMAT (1hO,8A,9HOELTA X J,6X,9HDELTA Y J»6X,9H0ELTA Z J) B 117 

ENO 8 118- 


SURfiOUTlNE ALP3ET (TIME, ALPHAD,BETAD»ETAD,CVP, THETA, PSI , PHI ) C 1 

COMMON / AB/ 1 1 1 , IWV,TX(50) * VVX (50) ,VVY (50) ,VVZ (50) ,H (50) , ALT (50) , V' C 2 

1X0 ( 6u ) , 9YO (50 ) C 3 

OR=. 01 79532925 C 9 

RO=57. 29577951 C 5 

CALL FTLUP ( T l ME , V X 1 , l , I I I , T X , V V X ) C 6 

CALL PILUP (TIME, VY1,1, III, TX,VVY) C 7 

CALL FTLUP ( T 1 M£ , VZl , 1 , I I I ,TX, WZ) C 8 

CALL FILUP (T IM£,Z, 1 * 1 1 1 ,TX,B) C 9 

CALL PTLUP <2 , 0X0* 1 , I *V, ALT , VXO) C 10 

CALL FTLUP ( Z , VYO , 1 , I WV, ALT- VYO) C 11 

VZ0=0. C 12 

PSlR=PSl*uw C 13 

PrlIR = PrtI°Oh! c 19 


33 



APPENDIX B - Concluded 


T tit l AP=TnETA»uP 

C 

15 

Al=SlN(PSlK) 

c 

lb 

X2=CC5 ( HS 1 ) 

c 

17 

A3-5[N(PhIk) 

c 

1(4 

A4=C0S<PnlP> 

c 

19 

X5 = S 1 :m ( T Ht’T AK) 

c 

20 

x6=C05 < TmETAk) 

c 

21 

B 1 =A6'>X2°VA1 ♦ Ab^Xl »VY1-XS°VZ1 

c 

22 

B2 = ( »?»XS»Xj-Al »X4> »VAl* (X2'->A4*Al«XS 0 X3) »VY 1 ♦ X6*X3* VZ 1 

c 

23 

( a2»X5»X4*a 1**X 31 *VAl *• <Al«xb°A4-X2*X3) »VY1 ♦X6*X4*VZl 

c 

24 

b^=Ao«A<;‘»'/A0+X'3'»XV<»VY0-XS‘'VZ0 

c 

25 

85=<X2«K5»A3-Xl»X4)»VAUMA2 i 'A4-t'Ai<»XS*X3> «VY0»X6»X3*VZ0 

c 

26 

b6=(»2'*A5<*X4 + Al»AJ.)'-*VXi>MAl»x5 t> X4-x2 o A3)*VY0*Xb<*X4**VZ0 

c 

21 

UP=3I-P4 

c 

28 

VP = r)2-P5 

c 

29 

kP -6 3-io 

c 

30 

ALPhAsATAimZ 

c 

31 

£ 1 As 4 r A-l? ( (SOhT ( ) »UM) 

c 

32 

eETa=ArAN2(vp»uP) 

c 

33 

ALPrtaO = ALPnA*«u 

c 

34 

B£ f AO=‘itT A*n>9 

c 

35 

ETAO=£TA»PO 

c 

36 

V2=UP»UP» VP»VP*wP»WP 

c 

37 

CVP=SQPT(V2> 

c 

38 

PETUw^ 

c 

39 

ENO 

c 

40 
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LANGLEY LIBRARY SUBROUTINE MATIN V 

Language : FORTRAN 

Purpose : MATINV solves the matrix equation AX = B, where A is a square coeffi- 
cient matrix and B is a matrix of constant vectors. The solution to a set of simul- 
taneous equations, the matrix inverse, and the determinant may be obtained. If the - 
user does not want the inverse, use SIMEQ for savings in time arid storage. For the 
determinant only, use DETEV. 

Use : CALL MATINV(A,N,B,M, DETERM, IPIVOT, INDEX, NMAX,ISCALE) 

A A two-dimensional array of the coefficients. On return to the calling 

program, A“* is stored in A 

N The order of A, 1 S N = NMAX 

B A two-dimensional array of the constant vectors B. On return to the 

calling program, X is stored in B 

M The number of column vectors in B. The expression M = 0 signals 

that the subroutine is used solely for inversion; however, in the CALL 
statement an entry corresponding to B must still be present 

DETERM Gives the value of the determinant by the formula 
DET(A) = (l0 100 ) ISCALE (DETERM) 

IPIVOT A one -dimensional array of temporary storage used by the routine 

INDEX A two-dimensional array of temporary storage used by the routine 

NMAX The maximum order of A as stated in the DIMENSION statement of the 

calling program 

ISCALE A scale factor computed by the subroutine to keep the results of compu- 

tation within the floating-point word size of the computer 
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Restrictions : Arrays A, B, IPIVOT, and INDEX have variable dimensions in the subrou- 
tine. The maximum size of these arrays must be specified in a DIMENSION statement 
of the calling program as A( NMAX, NMAX), B(NMAX,M), IPIVOT(NMAX), and 
INDEX(NMAX,2). The original matrices A and B are destroyed. They must be saved 
by the user if there is further need for them. The determinant is set to zero for a 
singular matrix. 

Method : Jordan's method is used to reduce a matrix A to the identity matrix I through a 
succession of elementary transformations l n , £ n _j, . . l\. A = I. If these trans 

formations are simultaneously applied to I and to a matrix B of constant vectors, the 
results are A - * and X where AX = B. Each transformation is selected so that the 
largest element is used in the pivotal position. (See ref. (a).) 

Accuracy : Total pivotal strategy is used to minimize the rounding errors; however, the 
accuracy of the final results depends upon how well-conditioned the original matrix is. 

Reference : (a) Fox, L.: An Introduction to Numerical Linear Algebra. Oxford Univ. 
Press, 1965. 

Storage : 54 2g locations. 

Subroutine date : August 1, 1968. 
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LANGLEY LIBRARY SUBROUTINE FTLUP 

Language : FORTRAN 

Purpose : Computes y = F(x) from a table of values using first- or second-order 
interpolation. An option to give y a constant value for any x is also provided. 

Use : CALL FTLUP(X, Y, M, N, VARI, VARD) 

X The name of the independent variable x. 

Y The name of the dependent variable y = F(x). 

M The order of interpolation (an integer) 

M = 0 for y a constant. VARD(I) corresponds to VARI(I) for 

I = 1, 2, . . ., N. For M = 0 or N i 1, y = F(VARI(1)) for any value 
of x. The program extrapolates. 

M = 1 or 2. First or second order if VARI is strictly increasing (not 
equal). 

M = -1 or -2. First or second order if VARI is strictly decreasing (not 
equal). 

N The number of points in the table (an integer). 

VARI The name of a one -dimensional array which contains the N values of the 
independent variable. 

VARD The name of a one -dimensional array which contains the N values of the 
dependent variable. 

Restrictions : All the numbers must be floating point. The values of the independent 
variable x in the table must be strictly increasing or strictly decreasing. The fol- 
lowing arrays must be dimensioned by the calling program as indicated: VARI(N), 
VARD(N). 

Accuracy: A function of the order of interpolation used. 
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References: (a) Nielsen, Kaj L.: Methods in Numerical Analysis. The Macmillan Co., 
c.1956, pp. 87-91. 

(b) Milne, William Edmund: Numerical Calculus. Princeton Univ. Press, 
c.1949, pp. 69-73. 

Storage : 430g locations. 

Error condition : If the VARI values are not in order, the subroutine will print TABLE 
BELOW OUT OF ORDER FOR FTLUP AT POSITION xxx TABLE IS STORED IN 
LOCATION xxxxxx (absolute). It then prints the contents of VARI and VARD, and 
STOPS the program. 

Subroutine date: September 12, 1969. 
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BASIC DATA 


Computer Card Format 

A typical basic data computer card format is as follows: 

Column 

1 shows whether the data are for a landmark or fiducial point and, if 

the latter, which fiducial system is being used: 4 is used for land- 
marks; 1 for frame corner fiducial points; and 0 for film sprocket 
hole fiducial points 

2 to 5 vehicle or test identification numbers 


7 


onboard camera identification number 


8 

9 to 10 

11 

12 

13 to 14 

15 to 19 
20 to 40 
41 to 47 

48 to 50 


blank 

landmark number for landmark readings; 00 for fiducial point 
readings 

film reading system identification number 

blank 

1, 2, or 3; indicates which of the required three readings (for aver- 
aging purposes) of each fiducial point or landmark is on the card 

accumulated point number, including fiducials 

blank 

x- value (in inches x 10" 3) of the landmark or fiducial point on the 
frame, measured horizontally on the projected image, with origin 
at right edge of image or right sprocket holes of film 

blank 
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Column 

51 to 57 y- value (in inches x 10-3) 0 f the landmark or fiducial point on the 

frame, measured vertically on the projected image, with origin at 
bottom edge of projected image or bottom sprocket holes of film 

58 to 63 blank 

64 to 67 frame number based on a convenient zero reference frame 

Only the data in columns 9 to 10 and after column 40 serve as inputs to the com- 
puter program. The remaining data are used for identification and gross reading or 
system error analyses only and may be altered at the user's discretion. For the system 
employed at Langley Research Center, the data through column 11 are input manually 
on dials located on the digitizer. The remaining data are input automatically and dis- 
played on the digitizer. All data are recorded by the typer as well as on the punched 
computer cards. 


Data Reading 

The following discussion is related to the procedures and system currently 
employed at the NASA Langley Research Center utilizing the Gerber film reading system 
and focuses on getting the basic data into a form usable by the computer. 

First a frame reference time must be established and the landmarks and features 
to be used must be identified. Next, sketches showing the position of the landmarks (and 
features) in relation to other distinct features (for example, mountain ranges, rivers, 
surface discolorations) on the film should be prepared. These sketches are useful in 
determining the general landmark location on the projected film image. 

A typical sequence in reading a frame of film, after the film has been installed in 
the film reader, is enumerated. These instructions are for the Gerber film reading 
system at the Langley Research Center and may be modified for other systems. 

(1) Set the frame counter to the proper value 

(2) Establish the image space axes system (x, y) origin (for example, frame corner) 

(3) Focus the frame (do not change the focus until the next frame) 

(4) Assure that the x and y values increase in the desired direction by moving 

the crosshairs 

(5) Aline the crosshairs on the first fiducial point (for example, frame corner, 

film sprocket hole) 
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APPENDIX E - Concluded 

(6) Assure that all dials on the digitizer are properly set and reset the x and y 

values on the digitizer to zero 

(7) Assure the readiness of the typer and card punch 

(8) Punch the fiducials (three readings for each of the four fiducials) 

(9) Adjust the digitizer dials as necessary for landmark readings and punch the 

landmarks (three readings each) 

(10) Check the printing from the typer for general correctness 

(11) Advance to the next frame to be read and repeat these steps. 

Finally, the cards should be interpreted and listed. A review of this listing should 
reveal any errors or blank cards which should be removed prior to processing. 
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APPENDIX F 


CONTROL DATA AND BASIC DATA INPUT DESCRIPTION 


Control Data 


For user's convenience, program VOPD has been divided into distinctive input and 
output groups. Following is a detailed description in computer terminology of the input 
for each group of control data. 

Column 

Parameter 

Description 

Units 

Input Group 1 

Format (2F10.3) 

25 cards 


1 to 10 

XRR 

Radial distance from frame 
center 

mm 

11 to 20 

RRR 

Lens radial distortion 

pm 

Input Group 2 

Format (80H) 

1 card 


1 

1 

Punch a 1 in column 1 


2 to 80 

TITLE 

Descriptive title, printed at top 
of each output page 


Input Group 3 

Format (215, 6F5.1) 

1 card 


1 to 5 

HI 

Number of data values in Input 
Group 7 


6 to 10 

IWV 

Number of data values in Input 
Group 4 


11 to 15 

SIM 

Multiplying factor to convert 
camera yaw angle (aj) to 



vehicle yaw angle {xp) 

16 to 20 S1A Addition factor to convert deg 

camera yaw angle (a^) to 
vehicle yaw angle (xp) 

PSI = SIM * SIG1 + SI A, 

SIM = +1 or -1 
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APPENDIX F - Continued 


Column 

Parameter 

Description 

Units 

21 to 25 

S2M 

Multiplying factor to convert 
camera pitch angle ( 'a 2 ) to 
vehicle pitch angle ( 0 ) 


26 to 30 

S2A 

Addition factor to convert 
camera pitch angle ^ 2 ) to 
vehicle pitch angle ( 0 ) 
THETA = S2M * SIG2 + S2A, 
S2M = +1 or -1 

deg 

31 to 35 

S3M 

Multiplying factor to convert 
camera roll angle ( 03 ) to 
vehicle roll angle (<p) 


36 to 40 

S3A 

Addition factor to convert 
camera roll angle (a 3 ) to 
vehicle roll angle (cp) 

PHI = S3M * SIG3 + S3 A, 
S3M = +1 or -1 

deg 

Input Group 4 

Format (3F10.1) 

IWV cards 


1 to 10 

ALT 

Altitude table, mean sea level 

ft 

11 to 20 

VXO 

North-South wind velocity 
table, + from South 

fps 

21 to 30 

VYO 

East -West wind velocity 
table, + from West 

fps 

Input Group 5 

Format (1615) 

(NFR + 3)/16 cards 


1 to 5 

NT 

Number of landmarks read on 
this film, also number of 
Input Group 9 data values 




APPENDIX F - Continued 


Column 

Parameter 

Description 

Units 

6 to 10 

NFR 

Number of frames read from 
this film 


11 to 15 

NT FR 

Number of Input Group 6 data 
values 


16 to 20 (etc.) 

NPT(I) 

Number of landmarks in each 
frame, etc., I = 1 to NFR 


Input Group 6 

Format (115, IF 10.1) 

NTFR cards 


1 to 5 

IFRM 

Frame number, does not have 
to correspond to frames 
read, but must encompass 
those of the basic data deck 


6 to 15 

TIMX 

Flight time of this frame 

sec 

Input Group 7 

Format (8 F 10.1) 

HI cards 


1 to 10 

TX 

Time table 

sec 

11 to 20 

XXX 

X-location of the vehicle 
from origin, + East 

ft 

21 to 30 

YYY 

Y-location of the vehicle 
from origin, + North 

ft 

31 to 40 

ZZZ 

Z -location of the vehicle 
from origin, + zenith 

ft 

41 to 50 

wx 

Vehicle velocity, + toward 
North 

fps 

51 to 60 

VVY 

Vehicle velocity, + toward 
East 

fps 
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APPENDIX F - Continued 

Column Parameter Description Units 


61 to 70 

VVZ 

Vehicle velocity, + toward 
Earth 

fps 

71 to 80 

H 

Altitude of vehicle, mean 
sea level 

ft 

Input Group 8 

Format (7F10.0) 

1 card 


1 to 10 

CC 

Camera focal length 

mm 

11 to 20 

PHI 

Latitude of the origin of 
Group 7 coordinate 
system, geodetic 

deg 

21 to 30 

THZ 

Longitude of the origin 
of Group 7 coordinate 
system, geodetic 

deg 

31 to 40 

RZ 

Distance from center of 
Earth to origin of Group 7 
coordinate system 

10 3 

41 to 50 

SIG1 

Initial estimate of camera 
Euler yaw angle 

deg 

51 to 60 

SIG2 

Initial estimate of camera 
Euler pitch angle 

deg 

61 to 70 

SIG3 

Initial estimate of camera 
Euler roll angle 

deg 

Input Group 9 

Format (1110, 5F10.1) 

NT cards 


1 to 10 

ITGT 

Landmark ID number 


11 to 20 

THJ 

Longitude of this landmark, 

deg 


degree part only 
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APPENDIX F - Continued 


Column 

Parameter 

Description 

Units 

21 to 30 

THJMI 

Longitude of this landmark, 
minute part only 

min 

31 to 40 

PHU 

Geodetic latitude of this 
landmark, degree part 
only 

deg 

41 to 50 

PHI J MI 

Geodetic latitude of this 
landmark, minute part 
only 

min 

51 to 60 

RJ 

Landmark altitude, mean 
sea level 

ft 

Input Group 10 

Format (115) 

1 card 


1 to 5 

NOUT 

Number of data values in 
Input Group 11 


Input Group 11 

Format (715) 

NOUT cards 


1 to 5 

IFM 

Frame number as read from 
film 


6 to 10 (etc.) 

JOUT(I) 

Landmarks to be ignored for 



this frame only, I = 1 to 6 
for as many as required 


If JOUT(I) = 100 this entire frame is ignored; do not remove its data from the 
basic data deck. 
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APPENDIX F - Concluded 


Column 

Parameter 

Description 

Units 

Input Group 12 

Format (1615) 

NT/16 cards 


1 to 5 

NGO 

Input group number to which 



program returns after pro- 
cessing all the frames of 
this job = 13 to stop further 
calculations 

6 to 10 (etc.) IOUT Landmarks to be ignored in 

computations for all frames 

Basic Data 

The basic data deck as punched from the film reader is input after the control data 
deck. It consists of three cards for each of four fiducial frame corners plus three cards 
for each of the NPT landmarks for a total of 12 + 3 * NPT cards. Following is a list of 
the essential information from the basic data deck required for program operation. For 
an explanation of inputs for all columns of the basic data cards see appendix E. 


Column 

Description 

Units 

9 to 10 

Landmark number 


41 to 47 

x -value of landmark on the film 

in. 

51 to 57 

y -value of landmark on the film 

in. 

64 to 67 

Frame number 



Only one copy of this basic data deck is required. For additional calculations with 
these data only these input groups, including and following the one specified by NGO, 
must be inserted behind this basic data deck. 
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APPENDIX G 


PROGRAM VOPD OUTPUT DESCRIPTION 

The output of this program is listed one page to a film frame and includes the 
following groups of data: 


Parameter 


Description 


The title card is printed at the top of the page. 


Output Group 1 


NPT lines 


L-MARK 


Landmark number 


LONGITUDE 


Landmark longitude 


LAT-GOCN 


Landmark geocentric latitude 


RAD-ER-KFT 


Earth center to landmark distance 


LAT-GODT 


Landmark geodetic latitude 


ALTITUDE -FT 


Landmark altitude, mean sea level 


XLAM 


XMU 


XNU 


Direction cosines of image space coordinate system 
relative to object space coordinate system (see 
appendix A) 


Output Group 2 


1 line 


X-location of vehicle, same as Input Group 8 XXX 


Y-location of vehicle, same as Input Group 8 YYY 


Z -location of vehicle, same as Input Group 8 ZZZ 
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APPENDIX G - Continued 


Parameter 

Description 

Units 

XPC 

X-location of film frame center, average of fiducials 

in. 

YPC 

Y-location of film frame center, average of fiducials 

in. 

CC 

Camera focal length 

mm 

Output Group 3 

NPT lines 


XCO FINAL 

X-location of landmark in image space (on film) 

in. 

YCO FINAL 

Y-location of landmark in image space (on film) 

in. 

Output Group 4 

NPT lines 


X 

X-location of landmark in object space 

in. 

Y 

Y-location of landmark in object space 

in. 

RDIST 

Landmark location from frame center, image space 

in. 

DELTA R 

Radial distortion of landmark due to camera lens 

in. 

Output Group 5 

1 line 


FRAME 

Frame number 


TIME 

Flight time 

sec 

SIG1 

Camera Euler angle, yaw 

deg 

SIG2 

Camera Euler angle, pitch 

deg 

SIG3 

Camera Euler angle, roll 

deg 

XP 

X-location of frame center in object space 

in. 

YP 

Y-location of frame center in object space 

in. 


L-9461 
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APPENDIX G - Concluded 


Parameter 

Description 

Units 

PSI 

Vehicle Euler angle, yaw 

deg 

THETA 

Vehicle Euler angle, pitch 

deg 

PHI 

Vehicle Euler angle, roll 

deg 

ALPHA 

Vehicle pitch angle of attack 

deg 

BETA 

Vehicle yaw angle of attack 

deg 

ETA 

Vehicle total angle of attack 

deg 

Output Group 6 

NPT lines 


RX(I) 

Error in X-coordinate from frame center 

in. 

RY(I) 

Error in Y-coordinate from frame center 

in. 

L-MARK 

Landmark number 


Output Group 7 

1 line 



The standard deviation of the residuals is listed here. 


50 



Example input and output data for several sample cases are presented in this appendix. 
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STANDARD DEVIATION OF RESIDUALS .11718 



HARK LONG ITl/OE LAT-GOCN RAD £R KM LAT-GOOT ALTITUDE-FT XLAM XHU XNU 

65 IGNORED eY USERS DIRECTION 

66 -106.54500 32.66227 20935.698 32.83833 10000.000 -.2227328 -.5279743 -.8195323 

67 -106.40833 32.65313 20929.698 32.82917 4000.000 .0349096 -.5390263 -.8415652 


APPENDIX H - Continued 


rO 

O' 

< <0 
h- • 

LU O 
eg 


LU • 
03 O 

rg 


< O' 

I r- 
cl 

3 • 

< in 

I 


rg 


m d* 

x » 

CL * 


CO 








o 









< g 









»- S3 









UJ • 



o 






X 00 



o 






►— m 



o 






1 



o 



OC 

g* 

* 




O -4- 




rg 

o 




U • 



< 

g* 

rg 




CJ 



K 

g- 


g 



H 



— J 

o 

o 

O' 






LU 

• 

• 

-A 






Q 

I 

1 

t— « • 



in 






03 g 



rg 






CL O' 



O CO 









CL rg 









> rg 








u 

• 








rg 

in 






rg 


o 







r> 







s> 

>0 

a. r*~ 


o 




►- 

o 

in 

> • 

X 

• 

in 



VO 

rg 

-4 

ro 

X 


rg 



*■* 

—4 

CO 


< 


O O' 



Q 

rg 

o 


X 


a m 




• 

• 

r- 

1 S3 h* 


X rg 



OC 


4- 

s> 

J O sO 


• 






CL 



r- 






X • 









in 


VO 









-J 









< 

OD 






rg 


3 

r*> 






r- 


a 

• 






O' 



O' 

o 

o 


ro 

o 

ro • 

s> in 

VO 

Nl h- 

-i r- 

r* 


O' 

g- 


— o m 

Ui 

—i 

< s) 

S3 


S3 

in 

^ GO 

g CO 

ac 

m 

z s> 

s> 

> 

rg 

in 

03 rg 

^ om 


*• 

~ s> 

S3 


—4 

g* 


> r- rg 

u. 


U. O' 



• 

• 

O 

OC O' co 

o 


r- 

* 


rg 

* 

g 

co ao 



O O' 

rg 




rg s> 

o o 

z 

00 

u • 

• 




lD • 

• • 

o 

• 

v rg 

s> 




00 

1 

M 

o 






J) m 


H' 

>■ O' 






l 


< 

ro 









fO 






g 


> 

1 






u O' 


Ui 


CJ 

o 


r- 

g 

o ^ 


3 


-t a 

o 


r** 

g 

• 

rg -s 


S3 

<4 o 

r- 


O' 

g 

03 <4* 

— eg fn 

Q 

in 

2 o 

si 

X 

g- 

03 

O' 

nvg 

ac 

• 

w o 

O 


O' 



— o in 

< 


il r- 

S3 


• 

• 


X »n o 

Q 

x r* 

rn 

S3 


r- 

O' 

□j 

ac cu O' 

Z 

O' 

o ** 

00 




x r- 

g ro 

< 

00 

o • 

• 




00 

o o 

K 

rg 

X H 

rg 




K • 

O) • • 

CO 

1 

.-4 





o 

— i 1 


Ui 







< 


3 

X 


X 



Ui 

3 


u 

oc 


ac 



X 

◦ 


*■* 

< 


< 



< o 



X 

x s> 

h- 

X 

S3 

r- 

oc <n 

(/) 


Ui 

i s> 

S3 

» 

S3 

s> 

ll r- 

Ui 


> 

3 


-J 




ac 



59 



APPENDIX H - Continued 


eg 

4 
< ro 
K- • 
UJ -4 

CM 


O 
4 eg 
►- co 

u; • 

CD O 
eg 


co 

<D 

CO 

ed 

o 

CD 


a 

S 




CO 




o rg 









o 

in m 




4 O' 





z 

O' gj 




X O' 





X 

g> in 




a. o 






>t *■< 




-j • 






-4 >r 

CO cu 




< m 
« 






• • 










1 1 








CL 






u 




UJ 


m 




«o 




X 

Z5 

to g> 




g> 




< 

x 

m eg 




x • 




o 

X 

O' o 




a. ro 






in O' 
to ro 




h- 




z 


in in 








UJ 


• • 




r- 




* 


1 1 




< o 




M 






h- ro 




-4 

x 

O' o 




UJ • 





<r 

r- 0 s 

o 



X ® 





-j 

«-4 O 




> to 




x 

X 

4* O' 

o 



1 






^ 4 

o 


at 4 4 







rg ro 

O 4 


eg o 





a 


rg o 

O • 


< 4 eg 





at 


• • 



»- 4 r» 

ao 




<■ 

h— 

1 



-1 O o 

in 




X 

u. 




UJ • • 

eg 




o 

• 

o a 



Q 1 1 

*-* • 




u. 

UJ 

o o 

in 



uo to 





Q 

O tJ 

rg 



CL. O' 






• • 

O ao 



•-g 





»- 

o o 

O. rg 







4 

M 

w o 

>■ rg 






ro 

1 

h- 

o o 

• 






O' 

> 

-J 

O 4 

in 



eg 



in 

< 

4 

«-g 




ro 



g> 






o g> 

a. o- 



o 


h- 




hOm 

> • 

X 


• 

h- 

a 

ro r* 

in 


i/> eg -g 

ro 

aC 



3 

o 

to -« 

eg 


NH ^ CO 


< 



-J 

o 

ro O' 

U O' 


Q eg o 


X 



CO 

i 

to eg 

a ro 


• • 

r* 

1 

g> e- 



h- 

CO co 

X eg 


at fO 4 


-J 

g> 4 



4 

• • 

• 



Q. -u 




O 

.j 

eg eg 

r- 



X • 




z 


fO fO 




in 



i/) 

•— « 









-j 

X 

x 








< 


* 

co co 

ao 



o 




> 


O' O' 

to 



ro 



Q 


QC 

4> 4 

• 



to 





UJ 

• • 

O' 

o o 

<o o 



O' O' 

to 



z m o' 

■m h- 

-j r- r* 

O' 4 


m 

UJ 


a O m eg 

—4 

< g> g> 

g» in 

rs ® 
/> ^ 

t —4 

rg ro 

at 


4 

m (j* O' 

ro 

z g> g? 

> eg in 

-mm 

in r— 



QC 

k- o o 

4 

« 4 4 

"4 4 


> 

eg ao 

u. 



U eg eg 

H 

u. o* ^ 

• • 

r- 

ct 

4 eg 

CJ 



UJ 


r- 4 

eg ^ 

3? 


ao ® 



z 

o 

OC ® ro 
•* rg ^ 

4 

00 

O O' eg 

u • • 



o o 

• • 

z 

a 


o o f- <o 

• 

>«N *> 




1 




m m 

g> 



ro 



K 


i 

</) g> gj 

> O' 



1 



4 


>- 

oc • • 

to 








< 

uj rg eg 

to 



_ ® 



> 


j uo cn cn 

1 



£ ** 



UJ 



J 


o o 

o- 4 




Q 





-J o o 

e- 4 

• 


rg ao 



UJ 

> o ro 

g> 

< o 

O' 4 


o o 

a 


Q 

flu to ro 

m 

z o g> 

X 4 ® 

O' 


q ao 

CL 



m oo 

• 

h o g) 

O' *-■ 


w 

4 

< 


> 

Q 4 o 


u. r- g> 

• • 


x m f 

a 



iu in g 

X 

ro g? 

r- ce 

UJ ^ 

a 

in gj 

z 


o a • • 

O' 

o <-* ® 


X h- 


4 ro 

< 


z a g> g> 

00 

o • • 


« ao 


O O 



o z o o 

eg 

X ^ eg 


H • 

CO 

• • 

l/) 


-j 

O *-4 

1 



O 

-j 

1 




~ 1 1 

UJ 



r- 

< 







X 

X 

UJ 





ac 


o 


a: 

X 

o 




4 

m g> 


< 

< 

< y 

M 




X 

g? gj gj 

X 

x g> e- 

x g? r*- 

ac rn 

00 




1 


UJ 

1 4> 'O 

i 4) g> 

u. r- 

UJ 




-J 


> 

-J 

-* 


aC 




60 



APPENDIX H - Continued 


eg 

*h 
< O' 
K- • 
LU -• 


< O' 

Uj • 
(S O 

in 


( 1 ) 

CO 

o$ 

o 


<13 


Cti 

CO 


m tn 
o © *g 

z © © 

x eg «*■ 

o' 

<0 N 
O' O' 

• • 
I I 

O' rg 
3 O' 4" 

x >o —* 

x © eg 

* ^ 
© r- 

't rg O 

I • • 

> I 

< 


< 0t 

X *+ 
% rg 


* 
in 
M O' 

X • 

a. 

© 


O' 

< <n 

h- ff> 



X 


r- 

O' 





UJ • 






< 


© 

-4* 

O 




X © 





<1 

_J 


© 


O 




h- © 





QC 

X 


© 

r- 

o 




1 





IU 



— « 

r- 

o 


oc r- 

o 






X 



pH 

© 

O 't 


<4 

«4* 






< 



o 

H 

u • 


< O' 

in 






u 



• 

• 

© 


h- *4* 


m 






h- 


t 

1 

pH 


-J O 

o 

m 






a. 






UJ • 

• 






QC 

« 


o 

o 



O 1 

1 

*4 « 





< 

UJ 


o 

o 

in 




m in 





© 

a 


o 

o 

rg 




a rg 








• 

• 

L_> m 




rg 






>- 


© 

o 

a 









a 

•— * 


u 

o 

> rg 








© 

< 

h* 


© 

m 

• 








r- 

o 



't 

M* 

in 




cn 




^g 

-j 

< 








4- 




rg 










a. r- 




O 


H- 






h- ® 

4- 

>• • 

X. 



• 

i- 

O 


o 

m 

m 


cO -< 

o 

cn 

QC 




Q 

o 


o 

<n 

fO 


M © 



< 




J 

o 


in 

m 

o © 


q <n 

in 


X 




© 

1 


e'- 

© 

a. o' 


• 

• 

m 

1 

© 

CM 





en 

ni 

X -« 


© m 

pH 

cn 

-J rg 

© 



< 


• 

• 

• 




cl -g 





O 

-J 


<n 

<n 

r- 




x • 





z 



<n 

m 





«n 




uo 














-J 


X 












< 


* 


© 

© 

r- 




^4- 




o 

> 



O' 

O' 

4- 




in 




Q 


QC 


O' 

O' 

• 




■n o' 






UJ 


• 

• 

m 

o o 

*• 

* 

o • 


in 

cn 

CO 



Z O' 

O' 

»M O' 

«j o r- 

* 

r* 




pH 

© 


O O rg 

CM 

in 

< o © 

O' 

© 

<n ® 

M 

»g 

O' 

QC 


< 

M4 

O' 

O' 

o 

2 0© 

>■ r- 

»h 

H 


O' 

rg 



QC 

1- 

o 

O 

rg 

*-• o © 

rg 



> 

in 

in 

U. 



O eg 

eg 

»«g 

U. f- Ul 

• 

• 

O' 

QC 

r- 

o 

CJ 



UJ 




O' © 

© 

rg 

CM ^ 


rg 

<n 



2 QC © 


r- 

O h- cn 



j> m 


O 

o 

z 


<-> .** «4* 


*h 

o • • 



-g * 


• 

• 

O 


oar- 

© 

• 

> © cn 



73 'O 



• 



o 


O’ 

© 

r- 




© 




h- 


i 

l/J 

*H 

o 

>• <n 




| 






»- 

QC 

• 

• 

© 








►4 


< 

LU 

cn 

m 

*■ 








> 


-j 

00 

m 

rf\ 

i- 




o ^ 




UJ 



J 




o o 

O' 

cn 





Q 







—j o r- 

© 

r- 

• 


O' 

in 



UJ 

> 

r^ 

O 


< o © 

o 

© 

in 


o 

m 

O 


o 

© 

© 

O 

rg 

z o © 

X ** 


rg 


© 

© 

CL 


3 


o 

in 

• 

•-H o © 


O' 

rg 

— • 

r» 

o 

< 


H- 

a 

<n 

o 

#-« 

ll *h 

• 

• 


X 

© 

o 

a 


P-4 

UJ 

rg 

m 

X «' 

© O' 

r- 

in 

uj r«- 

© 

© 

pH 

z 


O CL 

• 

• 

•g 

O rg rg 



X O' 


o 


< 


zoo 

© 


O • • 



M © 


o 

o 

>- 


a 

z o 

o 

m 

X o © 



h- • 

l/> 

• 

• 

LO 


-J o 


•H 


*-g 



| 

-J 


1 




i — • 

1 

1 

UJ 





< 





X 




-J 

X 

X 


UJ 

o 





ac 




o 

QC 

QC 


X 

a 





< 

O' 

© 

rg 


< 

< 


« o 

►4 





X 

© 

rg 

© 

X 

X © rg 

X © 

rg 

QC © 






1 




© 

1 rg © 

1 rg 

© 

X 

© 









> 


-j 



OC 





61 



APPENDIX H - Concluded 


o 
rn 
< o 

H- • 
UJ -4 


4 > 

< cn 

K CD 
Ui • 
CD CD 

in 



z> 

o 

4* 

4- 

in 





< 

A 





z 

4> 

O' 





X 

r- 





>< 

CM 

4* 





Q. 

4 






4* 

CJ' 





-J 

• 






<3 

r- 





< 

OJ 






O' 

o' 






•-4 






• 

• 

* 






4 






* 













OJ 

Oi 






o 






O' 

a 





M 

4* 





z 

ao 

O' 





X 

• 





X 

4> 

H 





a. 

^4 






4" 

m 






GO 






>o 

r- 






^4 




4* 


OJ 

o 










1 


• 

• 






4) 




> 



» 





< 

ao 




«x 








K 

m 





z 

CJ 

s0 





uj 

• 





< 

'D 

H 

o 




X 

4) 




< 

-J 

CJ 

O 

o 





00 




It 

X 

«o 

r- 

o 





1 




Ui 


H 

o- 

o 


0£ O 

H 






x 



00 

u >t 


4* 

in 






< 


o 

•4 

o • 


< O' 

m 






o 


• 

• 

o 


K- >T 

-4 


in 





t- 

1 

1 

~4 


— J o 

o 


4> 





u. 





UJ • 

• 


4- 




QC 

1 

o 

o 



o 1 

1 


• 




< 

Ui 

o 

o 

<M 




U> 

in 




03 

o 

o 

o 

4* 




CL 

Oi 






• 

• 

»J OI 





OI 





H- 

o 

o 

a. 4- 









a 

*-* 

o 

CJ 

> oi 









<. 

h 

m 

m 

• 









o 

-1 

* 

4- 

in 





OJ 




-j 

< 








4- 










o 

o 

a 

h* 





»- 





h- -* 

CD 

> 

• 

X 



►- 

O 

o 

rf> 

o 


OO O 

O 


r0 

a£ 



Q 

O 

o 

cq 

o 


►4 %0 

in 



< 



-J 

o 

in 

m 

O o- 


a *n 

in 



Z 



CD 

1 


43 

Q. 00 


• 

• 


O 

1 

oo 

OJ 


►- 


OI 

X -4 


ac m 

*4 



-J OI 

4> 


< 

• 

• 

• 




a. 

*4 




O 

-J 

m 

m 

r- 




X 

• 




Z 


ff> 

m 






in 




*-* 













X 

z 













X 

00 

ao 

O' 





o- 




> 


O 

O' 

O' 





o 





ac 

O' 

O' 

• 




rO 

't 





ui 

• 

• 

4* 

o o 

OJ 

4 

O 

• 


4> 

4“ 


z 

o* 

O' 

oj O' 

-j m o 

CO 

r- 


p4 

— • 

«i 

4) 


Q O 

04 

04 

in 

< m o 

*Q 

Oi 

i/> 

CO 

M 

m 

OJ 


< *4 

O' 

IT* 

CJ 

z m o 

>• oo 

O 

^4 


4" 

4 


0C h- 

o 

O 

04 

« ^ o 

oi 

4* 



> 

44 

4* 


O 

Oi 

OJ 

f* 

U. 0- o 

• 

• 


4> 

QC 

UI 

O' 


Ui 




o o 

o 

OJ 

rvj 

CO 


m 

m 


z a: 

ao 

r4 

CD 

O ao cn 



lD 

<n 


o 

o 


u *4 

* 

"4 

OI 

u • • 



*4 

• 


• 

• 


O Q 

r*» 

•o 

• 

> cd m 



A 

4> 



i 


O 

O' 

ao 

4- 





ao 





I i/> 

»4 

O 

> cry 





1 





H* QC 

• 

• 

4> 










< UJ 

m 

rn 

4* 





in 





-J 00 

m 

f0 





O 

>o 









o o 

r*» 

ro 


■r 









J O Q 

4 

ao 

n 

• 


ao 

OJ 


UJ > 

r* 

O 

ui 

< o o 

■4* 

^4 

UI 


4* 

O' 


Q CD 

D 

O 

o 

zoo 

X Oi 



OJ 

*-* o 

r4 


3 

o 

in 

9 

M o o 

m 



OI 

w 

m 

4* 


h- a 

m 

o 

PO 

U. HN 

• 

• 



X 

m 

O 


»— * UJ 

04 

m 

x «r 

o CD 

0- 

in 

Ui 

OJ 

QC 

4 

4* 


o a 

• 

• 

*4 

O 04 Oi 



z 

*■4 


^4 

4 


z o 

>0 

WO 

<■4 

U • • 



mm 

4) 


o 

o 


u z 

o 

O 

ro 

X O CO 



h- 

• 

<✓> 

• 

• 


J o 

H 

H 


<•4 




1 

-J 


1 


*■4 

1 

1 

UJ 






< 




X 



-J 

X 

X 


Ui 






aC 



CJ 

0 C 

a: 


Z 


a 




< O' 

00 

OI 

M 

< 

< 


< 

in 





Z sO 

OI 

*0 

X 

Z CO OI 

Z ao 

Oi 

CC 

4> 

c/) 




1 



UJ 

1 OJ NO 

1 Oi 

>0 

u. 


UJ 




-i 



J* 

-J 

-J 




CC 




62 


STANDARD DEVIATION OF RESIDUALS .02806 



REFERENCES 


1. Bendura, Richard J. ; Henning, Allen B.; and Smith, Richard E., Jr.: Vehicle Attitude 

Determination With a Single Onboard Camera. NASA TN D-4359, 1968. 

2. Brown, Duane C.: A Treatment of Analytical Photogrammetry. AFMTC-TR-57-22, 

U.S. Air Force, Aug. 20, 1957. (Available from DDC as AD 124144.) 

3. Hogge, John E.: Three Ballistic Camera Data Reduction Methods Applicable to 

Reentry Experiments. NASA TN D -4260, 1967. 

4. Lu ndstrom, Reginald R. ; Raper, James L.; Bendura, Richard J. ; and Shields, 

E. William: Flight Tests of Viking Parachute System in Three Mach Number 

Regimes: I - Vehicle Description, Test Operations, and Performance. NASA 
TND-7692, 1974. 

5. Bendura, Richard J. ; Lundstrom, Reginald R. ; Renfroe, Philip G.; and LeCroy, 

Stewart R.: Flight Tests of Viking Parachute System in Three Mach Number 

Regimes, n. Parachute Test Results. NASA TN D-7734, 1974. 

6. Whitlock, Charles H. ; and Bendura, Richard J.: Dynamic Stability of a 4.6-Meter- 

Diameter 120° Conical Spacecraft at Mach Numbers From 0.78 to 0.48 in a Simu- 
lated Martian Environment. NASA TN D-4558, 1968. 


NASA -Langley, 1974 L-9461 


63 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
WASHINGTON, D.C. 20546 


OFFICIAL BUSINESS 

PENALTY FOR PRIVATE USE $300 SPECIAL FOURTH-CLASS RATE 

BOOK 


POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 
451 



POSTMASTER : 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 


"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof." 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. Also includes conference 
proceedings with either limited or unlimited 
distribution. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include final reports of major 
projects, monographs, data compilations, 
handbooks, sourcebooks, and special 
bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 





